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ABSTBACl 


This document describes a computer program known as 
POTGRH which has been developed as an independent segment of 
the NASA Ames Three-Dimensional Linearized, Potential Flow 
Analysis System (POTFAN) and which is used to generate a 
panel point description of arbitrary, three-dimensional 
bodies from convenient engineering descriptions consisting 
of equations and/or tables. Due to the independent, modular 
nature of the program, it may be used to generate corner 
points for other computer programs. 
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1 IMTflOpUCTIQN 


This document describes Version 1 of a geometry 
generating code (POTGEH) which computes panel corner points 
and other geometrical data necessary to run the remaining 
segments of POTFAN, which is a program system for analyzing 
three-diaiensional, subcritical potential flows about 

arbitrary configurations. An overview of the POTFAN system 
is given by Medan (1976). In addition to being the first 
segment of POTFAN, POTGEM can be used to generate 

geometrical corner point data for other computer programs- 

Some of the important characteristics of the POTGEM 
program are the following: 

1. Complicated compoments can be handled with a minimum of 

input. Components that can be handled include 
complicated aircraft fuselages; thick wings with 
variable section, twist, and dihedral; wings with 
control surfaces; wind tunnels; and fan wakes. One or 
more components considered together constitute a 

configuration. 

2. In addition to computing panel corner points, the 
program computes position vectors to control points, 
vortex force sensing locations, and other useful data. 

3. Simple panel distributions are easily handled, yet the 
program does allow complex and non-uniform panel and 
control point distributions. 

4. The program allows components to be decomposed into 
segments. The geometrical description of the component 
and the independent variables may be different in each 
segment. 

5. The program has a general rotate, shift, and scale 
capability. 

6. The program runs under a command format which makes the 
program flexible, easy to use, and easy to modify even 
though the program is quite complex and versatile. 

7. Variable dimensioning is used so that oddly-sized 
problems can be handled without redimensioning. 
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8. Hachxne dependent language teatures have been generally 

avoided to make conversions to other computers 

relatively easy. 

9. The program has been liberally documented internally 
with comment cards to make it easy to modify. 

10. The program checks tor user input errors in many 

places. This makes it somewhat difficult to improperly 
run the program. 

11. The program is coded in FOfiTFAN IV. 

12. To date the program has been run on a UNIVAC 1108 

computer with 65k words of memory. 

13. The program can be executed in either batch or 

conversational modes. 

14. The program is based on a generalized coordinate system 
that reduces to Cartesian, polar, or spherical in 
special cases, yet is more flexible than any of these. 

One major feature lacking in the program is that 
because each component of a configuration is considered 
independently, the program cannot automatically put panel 
edges along lines of intersection with other components. It 
is necessary for the user to supply these intersection lines 
in the form of the VL (S) , VU(S), SL(V), or SU (V) curves that 
are defined in Section 3.2.1. It should be noted, however, 
that the program produces output (S and V at corner and 
boundary condition points as described in Section 6. 1) that 
would allow another program to be developed which could 
automatically calculate the intersection lines. This would 
allow POTGEM to then be rerun with the correct intersection 
lines. Therefore, part of the intersection problem has 
already been solved in POTGEM. It is expected that such a 
program will be developed in the future. 

Another disadvantage of POTGEM is that it cannot 
automatically distribute panels in relationship to surface 
curvature (e.g. , dense panelling in regions of high 
curvature). Furthermore, POTGEM cannot be easily modified 
to do this automatically. However, the output from POTGEM 
(see Chapter 6) is suitable for driving another program 
which could redistribute panels based on surface curvature. 

Finally, the cross sectional data (Section 3. 1. 1) that 
is input to POTGEM is not in parametric form. If the 
program had initially been designed this way, it would have 
teen somewhat easier to use and, in some cases, multiple 
segments (Section 3.1.7) would not have been reguired. This 
deficiency is expected to be corrected in version 2 of 
POTGEM. 
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2 PROBL EM TA SK B£SC ^IPTIQ M 


This section describes the basic specifications that 
guided the development of this program and the basic 
mathematical problems confronting the authors at the 
beginning of the task. 

The task that the program was reguired to perforin is to 
produce a file called the geometry file and containing the 
data described in Section 6. It was required that this file 
be created from a convenient engineering description of as 
general a component as possible. Furthermore, it was 
required that the program be flexible, easy to use, easy to 
modify, well documented, and easy to convert to other 
computers. 

There are basically two mathematical problems 
associated with determining a panel corner point description 
of a component from an engineering description. The first 
is to devise a method which will give the position vector of 
any point on the surface given the two independent panelling 
variables, which are called S and V. The second problem is 
to divide an appropriate region of the s-V plane into 
quadrilaterals with the corner points of the quadrilaterals 
corresponding either to corner points of the panels or to 
control points of the panels. It is important to realize 
that these two problems are completely independent. 
Therefore, the computer program handles these phases 
separately and the method of solution of either of the 
problems is independent of the other. 
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3 METHOD OF SOL OTION 


This section describes the solutions of the two 
loathe mat ical problems posed in the previous section. 


3.1 SURFACE REPRESENTATION 

As shown in Figure 3.1-1, the surface is described in 
part by a set of cross sections, an arbitrarily curved axis 
(which is not necessarily perpendicular to the cross 

sections) , and the angular orientation of the cross 

sections. This data together with a method for 
interpolating between cross sections completely defines the 
surface. Each of these subjects is discussed in detail 
below. Following this there are explanations of how the 
pieces are fit together to mahe a worieiug algorithm and of 
the multiple segment capability. 


3.1.1 Cross Sections 

The set of cross sections defining the component 
consists of one or more members. Each member of the set may 
consist of an open or closed curve, but the curves are 
restricted to lie in a plane. 

The independent variable in the cross section is V and 
the dependent variable is V2. V may be either y* , z», or 0 
and V2 nay be either z*, y* or R, respectively, where R = 
SQRT(y***2 ♦ z»*»2), e = ATAN(z'/y'i# and where y' = 
YVyPSCAl(S) and z* = ZVZ^SCAL(S). Here Y« and Z* are 
defined by Figure 3.1-1 and YPSCAL and ZPSCAL are arbitrary 
scaling factors. These scaling factors are functions of the 
other independent variable, S. 

The choice of which pair of cross section variables 
that nay be used is restricted only by the requirement that 
V2(V) be a single valued curve. If none of the three 
choices yields a convenient, single valued curve, then the 
component must be broken into two or more segments such that 
for each segment V2(V) is single valued. 'The multiple 
segment capability will be discussed further in Section 
3,1.7. The choice of V and V2 cannot vary from cross 

section to cross section within the same segment, but may 
vary from segment to segment. 
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Once an appropriate set of cross section variables has 
been chosen, it is necessary to consider how the cross 
sections can be mathematically described. Each cross 
section may he described either by the coefficients of a 
series expansion (e.g. , V2 (V) =R (©) -AO*COS (6) ♦ 

AI^COS (0) +B1*SIN (0) ♦. .. ) or by a set of data points together 
with a specification of an interpolation method. In the 
latter case, the data points need not be (V,V2) pairs, but 
can be (y*#z*)» (z',y*)# ot (0,fi) pairs since the program 

can internally convert whatever is given to it into (V,V2) 
|:airs. The method used to describe any cross section is 
independent of the methods used to describe other cross 
sections (i.e., a table may be used for one cross section 
and a series expansion for another) . More details on the 
mathematical description of these curves are given in 
Section 3.1.4. 

The maximum number of cross sections allowed is 
governed by the variable MXD05 as explained in subroutine 
GEOM. 


3.1.2 Arbitrary Axis 

The axis equation has been chosen to be the following 
parametric form; 

XAXIS = XAXIS(S) 

YAXIS = YAXIS (S) 

ZAXIS = ZAXIS(S) 

The functions XAXIS (S) , YAXIS (S), and ZAXIS (S) can be 
described using the same methods used for the cross 
sections. These methods are described in Section 3.1.4. 
Each of the three functions can be defined independently. 

In view of the general form of the axis equation, the 
independent variable S can be identified with several 
physical quantities. If, for example, XAXIS (S) = s, then s 
is the value of X along the axis. S could also be Y or Z or 
the arc length. Furthermore, if XAXIS = YAXIS = ZAXIS = 0, 
S may even be an angle (see Section 7.5 for an example). 
Usually, however, S is the value of X along the axis. 


3.1.3 Orientation of Cross Sections 

It is well known (Euler's rigid body theorem) that the 
rotation of any rigid body can be effected by a single 
rotation about some axis. In the present case, the body is 
to be identified with the cross section and the rotation 
refers to the angle through which the Y, Z plane must be 
rotated to make it parallel to the Y', Z' plane. 

If the amount of rotation is denoted by PHI, and the 
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axis of rotation has the components EX, EY, and EZ, then the 
rotation can be defined by a 3-by-3 matrix, T, whose 
eguation is the following: 



where 





EX* = EX/SQRT (EX**2 ♦ EY**2 ♦ EZ**2) (3. 1. 3-2a) 

EY* = EY/SQRT (EX**2 + EY**2 + EZ*^2) (3. 1. 3-2b) 

EZ» = EZ/SQRT (EX + + 2 ♦ EY**2 ♦ EZ**2) (3. 1.3-2c) 


Each of the quantities PHI, EX, EY,and EZ is considered 
to be a function of S and can be described using the methods 
discussed in Section 3.1.4. 

If the matrix T is known, but PHI, EX, EY, and £Z are 
not known, then PHI, E)^, EY, and £Z must be calculated since 
the program works only with PHI, EX, EY, and EZ and not with 
the individual cciuponents of T. To calculate these 
quantities in the general case, one must first determine the 
eigenvalues and eigenvectors of T. One of these eigenvalues 
must equal one. The components of the eigenvector 
corresponding to the unit eigenvalue can be defined to be 
EX, EY, and EZ. PHI can then be determined by working 
backwards through equations 3. 1 . 3-2 and 3. 1.3-1. 

This cross section rotation is usually nonzero only in 
three types of problems. The first consists of a wing with 
a twist distribution; the second consists of a highly 
cambered fuselage whose cross sections perpendicular to the 
axis are much easier to obtain than those parallel to the 
Y-Z plane; and the third consists of an axisymmetric , but 
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otherwise arbitrary, body. In the first case the axis of 
rotation is the spanwise axis and PHI is the twist anyie. 
In the second case the axis of rotation is perpendicular to 
the plane in which the cambered body axis lies and PHI is 
the arctangent of the derivative of the camber. In the 
third case the axis degenerates to a point and PHI becomes a 
spherical polar angle (see Section 7.5 for an example). 


3.1.4 Representation of Functions 

In the previous three sections a number of functions 
were introduced, but no mention was made of how these 
various functions can be described to the computer program. 
This will now be explained. 

To begin with, c^ll of the various functions will be 
treated in the same way and will be referenced with the same 
variable names. Thverefore each of the functions requires an 
index to uniquely identify it. This index is the integer 
variable IC. Table 3. 1.4-1 lists the correspondence between 
IC and the various functions defining the surface. This 
table also shows in parenthesis the defaults of the various 
curves. 


As just mentioned, all of the functions described in 
the previous three sections will be referenced by the same 
names. In particular, VAR1 will stand for S or V (depending 
on the value of IC) , and VAR2 will stand for XAXIS, YAXIS 
ZAXIS, PHI, EX, EY, EZ, YPSCAI, 2PSCAL, or V2 at any inpu 
cross section (depending on the value of IC) . 


In addition to the index, there is an integer function 
option, COPT, which identifies the basic type of description 
to be used. For example, COPT= 1 implies that the function 
will be determined by linear interpolation from an input 
table and COPT= -2 implies that the function will be 
determined by a power series expansion from an input set of 
coefficients. The above types of function definitions as 
well as others are implemented in subroutine XNTRP3 and 
associated subroutines. Table 3. 1.4-2 lists the 
correspondence between the function option number and the 
types of functions available in INTRP3. Note the following 
correspondence between IHTRP3 variables and the variables 
associated with the functions introduced in the previous 
three sections: 


lOPT 
XIN (♦) 
YIN (♦) 
PARAN (*) 
NIN 


COPT (IC) 
VAR1 (*,IC) 
VAR2 (*,IC) 
PARAM (♦,IC) 
NTAB(IC) 
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Therefore, in order to describe to the computer what 
any of the functions in the previous sections is, basically 
what one needs to do is to input the value of 1C and the 
variables COPT (IC) and NTAB(IC) and the arrays (VAR 1 (*, IC) ) , 
(VAR2 (*,IC) ) , and (PARAM (♦, IC) ) . Then the program will 
determine the function, when required, using subroutine 
INTRP3. 

The above variables are frequently all that one needs 
to consider when inputting any of the curves discussed in 
the previous section. However, in certain applications 
involving tables, some additional data manipulation 
capability is useful or required. This capability consists 
of three transformations that may be performed on the arrays 
(VAH1) and/or (VAR2) . The first transformation consists of 
the following general affine transformation: 

VAR1» = AFTRAN (1) ♦VAR1 + AFTRAN (3) ♦VAR2 + AFTRAN (5) 

VAR2* = AFTRAN (2) ♦VAR 1 + AFT RAN (4) *VAR2 ♦ AFTRAN (6) 

This has obvious usefulness for scaling and shifting tables. 

The second transformation is a more general 
transformation of the arrays (VAR1) and (VAH2) and includes 
an affine transformation as a general case. However, the 
most typical use of this second transformation is when the 
curve is a cross section curve (i.e. IC >11) and (VARl) and 
(VAR2) constitute a table. In this case the second 
transformation is used to transform (VARl) and (VAR2) so 
that (VARl) will be an array of V values and (VAH2) will be 
an array of V2 values. This second transformation would 
thus be very useful if, fcr example, it were necessary to 
use polar coordinates for V and V2 and the availaole data 
were in Cartesian coordinates. Both the first and second 
transformations are effected by subroutine TBAN2 as called 
by subroutines SRFINl and/or SRFIN2. 

The third, and final transformation is somewhat 

different than the first two. In the first place it only 
affects the array (VAR1). In the second place it has no 
effect on the meaning of the variables V and V2. That is 
(VARl) does not become an array of V values until after the 
first two transformations while the third transformation, 
although changing (VARl) , does not have an effect on the 
significance of V. Another way of stating this is if a 
certain set of program input data that did not specify the 
third transformation were used to describe a component, then 
the same input data only with the third transformation 
invoked would describe the same component. Thus the use of 
the third transformation does not require changes in the 
remaining input data. Also the third transformation will 
not affect the S and V values at corner points or boundary 
condition points. The only affect that the third 
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transformation has is on the quality of interfolated values, 
lo clarify hoi» this is possible, consider a two-dimensional 
blunt airfoil in the plane with the nose at the 

origin. The upper surface of this component can be 
described as 

Z»(Y*) = A*Y'**.5 ♦ B*y* + cn***1.5 ♦ ... 

Obviously an interpolation using polynominal spline fits 
would require many points near the nose to be accurate. Now 
let Y» = SQHT (Y*) . Then 

Z'(Y'*) = A ♦ B^Y« + CY'»**2 ♦ ... 

This curve does not require a large number of points tc 
result in an accurate interpolation. This is the function 
of the third transformation, namely to effect a 

pre-interpolation transformation of the independent variable 
to result in a much more accurate curve fit. 

The main situation in which this transformation would 
be used is on cross sections (i.e, V2 (V) curves) when such 
cross sections are blunt nosed airfoils. See section 7.6 
for an example. This transformation is effected by 

subroutine TRAN1, 

In summary, each of the curves required to define the 

surface can be input either as tables or coefficients and 

there are transformations available to manipulate tables and 
improve the accuracy of table interpolations. The input of 
these curves is accomplished with the SRII and SR12 commands 
discussed in Section 5.2. It should be noted that each of 
the curves may be input in a completely independent manner 
{e.g., PHI(S) may be described by a power series, YPSCAL by 
linear interpolation from a table containing five data 
pairs, V2(V) at SCS(1) by CODIM interpolation from a table 
containing ten data pairs, V2(V) at SCS(2) by a Fourier 
series, etc.). 


3.1.5 Interpolation Between Cross Sections 

A typical set of cross sections at which data is given 
is shown in Figure 3. 1.5-1. Each of the vertical lines in 
the S-V plane is a cross section and, therefore, by using 
the methods indicated in the previous section, each is a 
line upon which the dependent variable, V2, can be 
determined. Now the methods used to determine V2(3,V) in 
the remainder of the S-V plane will be described. 

Consider an arbitrary point, P, in the S-V plane and a 
horizontal line drawn through this point. The horizontal 
line intersects the given cross sections and at each 
intersection the value of V2 can be calculated. These V2 
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values together with the corresponding values o£ S 
constitute a table. Interpolation from this table is used 
to determine the value of V2 at the point P. 

This interpolation in the S-wise direction is done 
using subroutine INTBP3 and, therefore, all of the methods 
available in I11TRP3 are available for this interpolation. 
The method actually used is governed by the variable lOPTS 
and the array (PARAMS) that ate entered with the PAUL 
command (see Section 5.2). These variables are the same as 
lOPT and (PAHAM) in subroutine INTRP3, respectively. 

As is the case for the various functions described in 
Sections 3.1.1 - 3.1.3, there is a pre-interpolation 

transformation that can be invoked to increase the guality 
of the interpolation. As mentioned in the previous section, 
this transformation can be used with no other changes 
required to the input. This transformation is governed by 
STOPT and (PARST) that are entered with the PANL command. 
These variables correspond to lOPT and (PARAH) in subroutine 
TRAN1, Which performs the transformation. The most typical 
use of this transformation in this instance would be for a 
fuselage with both ends blunt and located at SI and S2. For 
this case, STOPT should be 4 and PARST (3) = (Si ♦ S2)/2 and 

PARST (4) = (S2 - Sl)/2, 

It should be noted that the interpolation between cross 
sections is done p rior to rotating the cross sections, prior 
to putting the cross sections on the arbitrary axis, and 
prio r to scaling the cross sections. This is consistent 
with”the way in which aeronautical structures are generally 
defined and, therefore, this results in an easy to use 
method. 


3,1.6 Summary of Geometry Algorithm 

In the previous sections various aspects of the method 
were explained. This section explains how these pieces are 
combined to make a working algorithm. 

There are two major phases. The first is the geometry 
input phase and the second is the actual calculation phase. 

In the geometry input phase the functions XAXIS(S), 
YAXIS (S) , ZAXIS(S) , PHI(S), EX (S) , EY (S) , EZ (S) , YPSCAL (S) , 
2PSCAL(S), and a set of cross section curves, V2(V), are 
defined according to the method described in Section 3, 1. 4 
and using the SRI1 and SRI2 commands. (Commands are 
described in Section 5.) Also the definitions of V and V2 
are established by .one of the commands POLS, CARY, or CARZ. 

Next comes the calculation phase in which values of X, 
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y, and Z on the body surface are detecaiined for given values 
of S and V. This phase is performed mainly in subroutine 
SURFAS. The determination of the actual values of s and V 
tor which calculations will be made is independent of the 
geometry definition method and is discussed in Section 3.2. 

Let (S,V) denote one of the given values of S and V. 
Then the first step is the calculation of the s-wise 
interpolation table (Section 3.1.5) for the given value of 
V. Next, the given S and values of S in the table are 
transformed according to the value of STOPT (Section 3.1.5) 
to make the interpolation more accurate. Then the 

interpolation is performed. As a result V2(S,V) is 

determined. Then y* and z* are calculated from ¥ and V2 
according to the definition of V and V2 (Section 3.1.1). 
This definition is stored in the variable VTYPE, which is 
established by the POLR, CARY, or CARZ command. Then the 
values of JtAXIS (S) , YAXIS(S), ZAXIS(S), PHI (S) , HX(S), 
EY(S), JiZ(S), YPSCAL(S), and ZPSCAL(S) 

are calculated. The variables y* and z' are multiplxed by 
YPSCAL(S) and ZPSCAL(S), respectively, to yield Y* AND Z' . 
The cross section rotation matrix, T, is calculated from 
equations 3. 1.3-1 and 3- 1,3-2 and, finally, X, Y, and Z are 
calculated from 

X(S,V) = XAXIS(S) ♦ T(1,2;S)*Y» (S,V) + T { 1 , 3; S) »Z • (S# V) 

Y(S,V) = YAXIS(S) + T(2,2;S)*Y* (S,V) + T (2 , 3; S) *Z » (S# V) 

Z(S,¥) = ZAXIS(S) ♦ T(3,2;S)=<^Y' (S,V) T (3, 3; S) *Z • (S, V) 

The last two terms in the above equations represent the 

rotation of the cross section to its final orientation. 

In addition to being used tc find X, Y, and Z values of 
corner points, the above method is used to determine the X, 
Y, and Z values of boundary condition points and also to 
calculate the unit normals. This is in contrast to many 
existing programs, which determine boundary condition points 
and unit normals from the corner points. The latter method 
IS generally not as accurate. 

The calculation of the unit normals will now be 
discussed. Consider a point P in the S,7 plane. Let the 
points A, B, C, and D be arranged around P in the manner 
shown in Figure 3. 1.6-1. The program calculates the 
position vectors (RA, RB, RC, and RD) to each of these 
points, calculates the cross product of RC-RA with RD-RB, 
normalizes the result, and calls it the unit normal. If 
UNEPSS and UNBPSV are small enough (but not too small), this 
method is generally more accurate than using the panel 
corner points. It should be noted that a consistent sign 
convention has been applied so that the unit normal will lie 
on the same side of the surface as a vector in the 
Nl-cross-N2 direction. If S and V are either both 
increasing or both decreasing functions of the Ml and N2 


3-8 



Indices, then the unit normals will be in the S-cross-V 
direction. 


3.1.7 Surface Segmentation 

A component may be divided into a number of segments. 
There are several reasons why this is done. 

In the first place the geometry may be naturally 
segmented. For example the NASA Ames l2.2m-by-24. 4m 
(40*x80'} wind tunnel has a cross section in the shape of a 
square with semi-circles on each side. (Figure 7.3-1 shows 
one side of this configuration.) The top and bottom 
surfaces can be easily and exactly described in a Cartesian 
coordinate system and one side can be easily and exactly 
described in a polar coordinate system with origin 
coinciding with the center of the corresponding serai-circle. 
(See Section 7.3 for an example.) 

A second reason for segmentation rests in the fact that 
the cross section curves, V2 (V) , must be single-valued. An 
example of where this requirement necessitates segmentation 
is a thick wing section. In this case two segments are 
required. One is the upper surface and the other is the 
lower. 

A third reason is that the component may be too 
complicated to handle as a single segment. That is, there 
may be too many cross sections and/or table entries to fit 
in the program at once. 

A final reason is that there may be certain lines on 
the component along which panel edges must be constrained to 
lie. An example of this is a wing planform with a crank. 
In modelling such wings it is best to have panel edges at 
the spanwise location of the crank. 

In addition to the above examples, some components may 
have to be segraented for wore than one reason. A thin wing 
with a deflected control surface, for example, may be easier 
to describe in segmented form and also panel edges should be 
made to lie along the hinge line. 

Despite the differing reasons for segmentation, the 
program handles all cases of segmentation in the same 
manner. Different aspects of this treatment are discussed 
in the following paragraphs. 

The total number of segments and the number of panels 
in each segment must be established once for each component 
before any panelling is done. This is accomplished by a 
DSEGMENTS command (Section 5.2). Note that this command is 
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soMewhat different from th^> geoaetcy definition coauandS/ 
the SEGMENT coanand, and the PANI comuand because it refers 
to the entire conponent, whereas the others refer to a 
segue nt (except the geoaetry definition coaDands need not be 
repeated if the geoaetry description remains the same). 

Prior to panelling each segment, the segment must be 
identified (see the SEGMENT command. Section 5.2), the 
boundaries of the segment in the S,V plane must be defined 
(see the SL, SO, VL, and VU commands. Sections 3.2 and 5.2) , 
and the distribution of corner points and boundary condition 
points along the boundaries (see the SLBC, SUBC, VLBC, and 
VUBC commands. Section 5.2) must be established. Also, 
prior to panelling each segment, those aspects of the 
geoaetry that are different from those of the previous 
segment must be redefined. 

Because of the fact that boundaries between segments 
may represent situations where the actual geometry is 
discontinuous (e.g., between a wing and the side edge of a 
control surface) the program inserts a pseudo row of panels 
in between each segment. In many cases it is desirable to 
eliminate these rows. This is accomplished by the NBI1 and 
NBI2 commands (Section 5,2). 


3.2 PANELLING 

In the previous section, a method was described for 
determining the position vector to the surface given the two 
independent panelling variables, S and V. This section will 
now describe how the values of S and V are determined. For 
each segment (Section 3.1.7), the locations of corner points 
are determined from the boundaries of the segment, the 
inter|Sections of the corner and boundary condition grid 
lines with the boundaries of the segment, and a method for 
constructing grid lines from the boundaries and intersection 
points. Each of these topics is discussed in detail in the 
following sections. 


3,2. 1 Segment Boundaries 

All segment boundaries in the S-V plane are in the form 
of four arbitrary functions SL(V), SU(V), VL(S), and VU(S) 
as shown in Figure 3.2. 1-1. Although in most applications 
the four curves are straight lines, there are situations 
where more general curves are desirable or essential (e.g. 
on a fuselage where a wing intersects). 

Each of these four functions can be described to the 
program in basically the same manner that the functions 



describing the geometry can (Section 3.1.4). That is, 
subroutine INTRP3, vith all o£ its various options, is used. 
These functions are input using the commands SL, SU, VL and 
VU (Section 5.2). 

As will be seen later, one of the necessary steps in 
determining S and V values in the interior of the segment is 
the calculation of the intersections of the VL and VU curves 
with the SL and SU curves. The method that the program uses 
for this calculation is an iterative one that cannot be 
guaranteed to converge for all cases. This iterative method 
is described in subroutine GRID. A sufficient condition for 
convergence at the intersection of the VL(S) and SL{V) 
curves is that 

ABS(d(VL)/dS) ♦ABS(d(SL)/dV) <1 

Similar conditions hold for the other three intersections. 
In the extreme cases for which the above conditions do not 
hold, the intersections must be determined a priori and 
input to the program with a PANL command. 

The notation for these intersections is defined in 
Figure 3-2. 3~1 (e.g., (SSUVL, VSUVL) is (S,V) at the 

intersection of SU (V) with VL(S)). 


3.2.2 Grid Line Intersections with Boundaries 

Grid lines are lines in the S-V plane along which 
corner points or boundary condition points can be iocated. 
nore specifically, corner points and boundary condition 
points are located at the intersections of grid lines 
extending nominally in the S direction with those extending 
nominally in the V direction. These grid lines are defined, 
in part, by their interesections with the boundaries of the 
segment. This section explains terminology related to these 
intersections and how the intersections are input. 

A nondime nsional system is used to describe the grid 
line and boundary intersections- This system is illustrated 
in Figure 3.2. 2-1 for the VL (S) curve. The intersection of 
VL(S) with SL(V) is, by definition, at XGP--1 and the 
intersection of VL(S) with SU (V) is at XGP=+1- Furthermore, 
XGP, by definition, varies linearly with S. In the program 
the intersections of the V-wise corner point grid lines with 
the VL(S) curve are contained in the array (XGPVLC) and the 
intersections of the V-wise boundary condition point grid 
lines with the VI (S) curve are contained in the array 
(XGPVLB) . Similar noration applies to the other grid line 
boundary intersections (e.g., (XGPSUB) contains the 
intersections of the SU (V) curve with the S-wise boundary 
condition point grid lines) . 
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The program does not assume that any of the elements of 
the XGP arrays are equal to -1 or ■♦•1 (i.e., corner points 

need not lie on the segment boundaries). Also the program 
allows, for example, XGPSLC (I) tc be unequal to XGPSUC (I) , 
although the program will assume that XGPSUC (X) is equal to 
XGPSLC (I) unless instructed otherwise. 


The intersection point arrays (i.e,, (XGPSLC). etc.) 
are defined for the program using the SLBC, SUBC, VLBC, and 
VUBC commands. Although the arrays may be input explicitly 
with these commands, the usual procedure is to select one of 
a set of predefined rules and have the program calculate the 
values. These details are governed by the variable lOPT 
that is input with the SLBC, SUBC, VLBC, and VUBC commands. 
If lOPT is zero, then the intersections are expected as 
input. otherwise, the program will calculate the 
intersections. The correspondence between some valid lOPT 
values and the intersections calculated is shown in Figure 
3. 2.2-2 for the case of four panels and lOPT greater than 
zero. The spacings for XOPT less than zero are reversed 
from those of positive lOPT. For example, I0PT=-4 gives 
finer spacing near XGP=+1. For any ncnzero lOPT the XGP 
arrays will always be in ascending order. The computation 
of corner point and boundary condition control point 
spacings is done in subroutine XPANCP. 


As will be seen in the next section it will be 
necessary to calculate the S and V values at the 

intersections of the grid lines with the boundaries. The 
notation for these intersections is shown in Figure 3.2. 3-1. 
The characters C and B standing for corner points and 
boundary condition points, respectively, have been dropped 
because the procedure is identical for both. The 

intersections of the S-wise grid lines with the SL(V) 
boundary curve are given by the following equations: 


VGPSL (K2) = ( 1-XGPSL (K2) ) ♦VSLVl/2 ♦ 
(UXGPSL (K2) ) ♦VStVU/2 


(3.2. 2-1) 


SGPSL(K2) = SL (VGPSL (K2) ) (3.2. 2-2) 

In the above equation SL (VGPSL (K2) ) is calculated using 
subroutine INTBP3 and the data entered with the SL command. 
Similar equations are used for the other three grid line and 
boundary intersections. 


3.2.3 Intersections of S-Mise and V-Hise Grid Lines 

As mentioned in the previous section the nondimensional 
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description of the intersections contained in the arrays 
(XGPSLC) , (XGPSLB) , etc,, only defines the grid lines in 
part. Namely, these arrays define the grid line 
intersections with the boundaries of the segment and they do 
so nondimensionally . This section completes the definition 
by explaining how grid lines are extended to the interior of 
the segment and how their intersections with other grid 
lines are calculated- 

The method used is identical for corner point grid 
lines and boundary condition point grid lines, and 
therefore, the characters C and E in the arrays describing 
the intersections will be dropped for the remainder of this 
section. The particular value of (S,V) to be determined is 
denoted by (S (K1 ,K2) ,V (K 1 ,K2) ) where K1 is an index that 
varies in the S direction and K2 is an index that varies in 
the V direction. All of the various symbols introduced to 
this point are shown in Figure 3. 2.3-1, 

A simple, yet effective way to determine S(K1,K2} would 
be to consider it a weighted average of the S values at the 
intersections of the s-wise grid line with the SL (V ) and 
SU(V) curves. Mathematically this can be stated as 

S (K1,K2)= (1-FS (K1,K2) ) ♦SGPSL (K2) / 2 ♦ 

(3.2. 3-1) 

(1 + FS (K1,K2) )*SGPSU(K2)/2. 

FS(K1,K2) is a nondimensional number that is equal to -1 
when the V-wise grid line coincides with SL(V) and is equal 
to +1 when the V-wise grid line coincides with the 3U(V) 
curve. A similar expression is used for V(K1,K2): 

V (K1,K2) = (1-FV (K1,K2) )*VGPVL(Kl)/2 + 

(3. 2. 3-2) 

(UFV (K1,K2) )*VGPVU (Kl)/2. 

This raises the question of how to determine F3(Kl,K2) 
and FV(K1,K2). The answer is shown in Figure 3. 2.3-2. This 
figure shows that FS (K1,K2) and FV (K1,K2) are at the 
intersection of two straight lines drawn from (XGPVL (K1 ) ,-1) 
to (XGPVU (K1) , +1) and from (-1, XGPSL (K2) ) to (+ 1 ,XGPSU (K2) ) . 
In other words FS(K1,K2) and FV(K1,K2) are the simultaneous 
solution to 

FS (FV (K1,K2)) = FV (FS (K1 ,K2) ) . 

More explicitly^ FS (K1,K2) and FV(K1,K2) are the simultaneous 
solution of 

FV(K1,K2) = (1-FS(K1,K2)) ♦XGPSL(K2) /2 ♦ 

(3,2. 3-3) 

(1+FS(K1,K2)) ♦XGPSU (K2) /2 
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and 


FS(K1,K2) = (1-FV(K1,K2))*XGPVL{K1)/2 ♦ 

(3. 2. 3-4) 

(1+FV (K1 ,K2) ) ♦XGPVU (K1) /2 

Due to their linearity, the above equations are easily 
solved analytically. Subroutine GRID contains the actual 
equations. This completes the description ct hov S and V 
values are determined. 

In summary, the V1(S), VU(S), SL(V) , and SU(V) segment 
boundary curves are defined by the VL, VU, SL, and SO 
commands, respectively. The grid line intersections with 
the boundary curves are defined by the VLBC, VU6C, SLBC, and 
SUBC commands. Then, when a PARL command is given, the 
program executes the following: (1) It calculates the 

boundary curve intersections ( {SSLVL, VSLVL) , etc.); (2) 
Using equations 3.2.2- 1 and 3. 2.2-2 and similar equations, 
it calculates S and V values on the boundaries; (3) It 
solves 3. 2.3-3 and 3.2. 3-4 simultaneously; (4) Finally it 
calculates S and V from equations 3. 2.3-1 and 3. 2. 3-2. The 
program executes the preceding (as required, i. e,, redundant 
calculations are avoided) for all of the corner points and 
all of the boundary conditions points and, also under the 
EARL command, determines the position vectors to the corner 
points and boundary condition points and the unit normals at 
the boundary condition points using the methods described in 
Section 3.1. 

An example of an actual network of grid lines generated 
by the above method is shown in Figure 3. 2. 3- 3. 
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a PBOGRAM DBSCBI PTIO N 

To a great extent the description of the inner workings 
c£ the program has been relegated to conment cards in the 
FOBTBAN source decks. This includes descriptions of the 
functions of the subroutines and their input and output. 
The remainder of the section presents relevant descriptive 
data which could not effectively be placed on comment cards. 


a. 1 CALLING STBUCTUBE 

Figure 4.1-1 shows the subroutine calling structure. 
Table 4.1-1 shows the calling structure in a different 
format. 


4.2 COMMON BLOCKS 

Table 4.2-1 shows the common blocks. These common 
blocks are the same length in every program in which they 
appear. 

4.3 LOGICAL UNITS 

Table 4.3-1 summarizes the logical units (tape, disks, 
or drums) which the program uses. 


4.4 MEMOay REQUIREMENTS 

Without the 25 arrays dimensioned in the main program 
and without using overlays, the POTGEH program reguires 
approximately 45,000 decimal words of core storage. This 
reguirement includes all system subroutines and internal 
symbol dictionaries and was determined on the INFONET Univac 
1108 operating system. The size of the 25 arrays must be 
added to this number to determine the total amount of 
storage reguired by the program. An overlay structure and 
compilation without internal symbol dictionaries can be used 
to decrease the storage reguired. 
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^.5 SYSTEM DEPENDENT SOBfiOUTlNES 


Subroutines FlbEND, OPENIi, and TIMES! all call system 
dependent subroutines. Therefore, they all generally must 
te modified when the program is used on a new computer 
system. 


4.6 BESTRXCTIONS AND LIMITATIONS 

The most important restriction regarding the POTGEM 
program is the limitation on the total number of corner 
points. This is limited by the dimensions of the 25 arrays 
in the main program. If these dimensions are not 
sufficient, then they must be increased, the main program 
must be recompiled, and the program must be relinked- 

There are also limitations on the maximum number of 
corner points in the N1 or S direction and in the N2 or V 
direction, the maximum number of cross sections, the maximum 
number of table entries, etc. These limitations are 
discussed in more detail in subroutine GEOH, instructions 
for changing the maximum limits are also given in GEOH. 

The program itself checks for violations of the above 
restrictions, so there is no a priori need for the user to 
worry about them. 
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5 OPERATING I NSTRUCTION S 

The purpose of this section is to provide the user with 
inforfflation necessary to execute the prograo. Instructions 
tor modifying the program ace also given. 


5.1 GENERAL DATA INPUT CONSIDERATIONS 

The program is designed to use commands as the basic 
form of input to control the program flow. These commands 
consist of four letters placed in the first four columns and 
are recognized as keys that cause the program to perform 
particular operations. These operations consist of reading 
input, writing output, or calculations, or a combination of 
all three. After the operations are performed, the program 
flow returns ,to the beginning of the program and reads the 
next command. This continues until a STOP command is 
encountered, whereby the program terminates. Any command 
input record whose first four columns are left blank is 
considered a "comment” command. In the conversational mode, 
any command that is not recognized by the program is prxnted 
and program flow is returned to the next command without any 
operations being performed. In the batch mode an 
unrecognizable command causes program termination unless the 
variable CONTIN has been made .TRUE. by a preceding DATA 
command. Following each command some data must usually be 
entered. This data is prescribed in either NAMELIST or 
regular formats. 

In the batch mode each command line that is read in is 
printed out before any action is taken on it. The entire 80 
columns are printed even though only the first four columns 
comprise the command. This allows the output to be 
documented with "comment" commands and helps to pinpoint 
sources of errors. 

Aside from commands, most of the input data is in 
NAMELIST format. The program has been coded to take maximum 
advantage of the way in which NAMELIST works. In 
particular, only the specific data that is actually required 
need;s to be input, the data elements can be in any order, 
and !(except where noted otherwise) data that is entered with 
one Command need not be reentered with succeeding commands 
unless it is to be changed. To effect this last advantage, 
intermediate arrays are used for some commands (SHIl, VLBC, 
VUBC, SLBC, SUBC, VL , VU, SL, and SU) . These intermediate 
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arrays allow sinilar curves stored in different core 
locations to be input with a BiniDum of data. 

Other input data is in regular format foc». Integers 
are always input under a 1615 format and floating point 
numbers under an 8F10.0 format. 


5.2 INPUT DESCRIPTION 

Detailed descriptions of the various available commands 
are given in the listing of subroutine GEOM. This 

subroutine is arranged in a number of sections. The first 
section consists of general commentary and a summary of the 
available commands. The next section is the specification 
section (sets data types, dimensions, common blocks, etc.) . 
The next section sets the initial default values for those 
variables having default values. The next section is the 
command read and branch section. By examining this section, 
one can determine to which point control is transferred for 
specific commands. Each of the following sections contains 
detailed commentary and code necessary to effect each 
command. Thus by determining the command transfer point 
from the command read and branch section, one is led to the 
place in the listing that contains the detailed description 
of the command and further input required for the command. 
Since the listing contains the detailed command and input 
descriptions, they will not be repeated here. 

Instead, the most typical ordering of commands and some 
general comments will be given: 

1 - TITLE 

2. DSEGMENTS 

3. Geometry definition — SRI1, SRI2, POLR, CARY, and 
CARZ. 

4. SEGMENT. Not required if N SEGS-NSEGV= 1 for the first 
segment to be panelled. 

5. Segment boundary definition SL, SU, Vi, and VU. All 

four commands are frequently not required for the 
second and subsequent segments. 

6. Corner and control point distribution definition -- 
SLBC and VLBC commands. Use Figure 3. 2. 2-2 as a guide 
to available distributions. SUBC and VUBC are not 
usually required. 

7. GRID. Prints distributions determined in step 6. 

8. PANL. 
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9. Repeat 3 thcough 8 as required for the remaining 

segments. The segments may be considered in any order. 
However, for each specific case, there is generally an 
optimal order that will minimize the amount of data 
required. Unless indicated in the detailed 

instructions, it is not necessary to repeat data that 
does not change. If all of the data input with any 
SRI1, SRI2, SL, SU, VI, or VU command remains 
unchanged, the command may be left out. 

10. NRI1 commands if there are multiple segments in 
V-direction and component is continuous. 

11. HRI2 commands if there are multiple segments in 
S-~direction and component is continuous. 

12. Singularity Definition — BCFLAG, DSFIAG, and UVW. 

13. ROSS. May be given more than once. 

14. FINISH. Enter variables in Section 6.1 that are marked 
with an asterisk (♦) . 

15. STORE 

16. PRINT 

17. STOP or compute a new geometry beginning with 
INITIALIZE followed by steps 1 through 16 or compute a 
modified geometry starting with an XINIT command 
followed by a subset of steps 1 through 16. 

It needs to be stressed that the above outline is not 
to be considered a rigid one. Ihe commands, except where 
noted, may be given in a different order. Also not all of 
the commands need to be given in all applications. 

Also it should be noted that segments do not need to be 
panelled in any specific order. 


5.3 SYSTEM COMTHOL CARDS 

This section describes the control cards that are 
necessary to run POTGEM on the various computer systems that 
have or are being used to run it. . 


5i3. 1 INFONET Uni vac 1108 System 

Since this system allows automatic file definition 
commands determined from the file identification number 
(Apppendix A) entered after the STORE command, the only 
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control card required is the name of the main program which 
is POTGEM/POTF. All geometry files created by the STORE 
command will show up in the user's LIBi library with names 
identical to the file identification number and a version 
identifier of GM- Thus, for example, inputting ID-1U23 
after the STORE command will create a file named 
1023.GH/LIB$. The content of these files is described in 
Section 6. 

In addition to POTGEM/POTF, the user may also want to 
switch IN$ and OOT$. 


5.4 PROGRAM MODIFICATIONS 

There are a few modifications that receivers of this 
program might typically want to change. These are described 
in the following sections. 


5.4.1 Additional Curve Fitting Capability 

If the curve fit routines provided are net adeguate, it 
is a simple matter to add new ones. All that is required is 
that the input can be made compatible with the input to 
subroutine INTRP3. If it can, then a new option can be 
inserted into INTRP3. An examination of the INTRP3 listing 
will make it obvious how this can be easily done. The only 
changes required to POTGEM are some simple additions of 
cards in INTRP3. 


5,4.2 Additional Types of Panel Distributions 

New panel and control point distributions can be easily 
added by simple modifications to subroutine XPANCP. This 
subroutine calculates the distributions demanded with any 
SLflC, VLBC, SUBC, or VDBC command. An examination of the 
listing of XPANCP will make it obvious how the modifications 
can be easily done. The only changes required are some 
simple additions of cards in XPANCP. 


5.4.3 Increasing Array Dimensions 
See Section 4.6. 

5.4.4 Use of POTGEM not in Conjunction with Other POTFAN 
Modules 

Persons wanting to use POTGEM to generate corner point 
descriptions of geometries for use by their own programs can 
do so without modifying POTGEM. However, by eliminating 
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uiineeded arrays and the code that conputes them, savings in 
CPU time and core memory can be realized. This section will 
explain how to eliminate the following arrays: 

(UVWX) , (UV«Y) , (UVWZ) , (BCFLAG) , (DSFLAG) , (SSFLAG) , 
(UNX) , (UNY) , (UNZ) , (DA) , (XSl) # (YSI) , (ZSl) , 

(XS2) , (YS2) , and (ZS2) . 

First of all these arrays should be eliminated in the 
main program. In the CALL GEOtt statement replace each of 
the array names with XCP (N) Where N is an integer larger 
than the core memory. The purpose of this is to ensure that 
a fatal execution error occurs in case the program should 
happen to try to compute one of these arrays. It is not 
necessary to remove the DIMENSIONS for these arrays in 
subroutine GEOM. 

The next step is to remove the code that calculates 
these arrays. An examination of a compiler cross reference 
map of subroutine GEOM will indicate which code can be 
removed. In particular note that the second CALL PANL2 
after a FANL command can be removed and so can the calls to 
FSENSI and FSENS2. Furthermore, note that some entire 
subroutines called by POTGEM can be removed. Also make 
appropriate modifications to (LOG) (see Section 6.2). 

The final step is to remove the code ttiat writes out 
these arrays. These are written out in subroutine STOKGM 
and PRMTGH. An examination of compiler cross reference 
listings of these subroutines will indicate clearly what to 
remove. Note that it is not necessary to eliminate the 
arrayi^ from the argument lists. 

If these changes are made, then like changes must be 
made in the auxiliary subroutines READGM and BD6HA and in 
the EDXTGM and PLOTGH programs. (see Section 8). 
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6 PROGRAM OUTP UT 

Output from POTGEM consists of line printei: output and 
geometry files. The line printer output is meant to be self 
explanatory and will not be discussed further. 

The geometry files are created according to the 
procedure in Appendix A and conform to the for it in 
Appendix B. Control cards for managing the geometi^ files 
are given in Section 5.3. Each geometry file consists of 8 
to 18 binary records (i.e., each is created by 8 to 18 write 
statements of the form WRITE (NIG) ...data...). The data 
contained on the geometry file is explained in detail in 
Section 6.1 and is summarized in Section 6.2. 


6. 1 Detailed Description of Geometry file 

All of the geometry file data on the most general 
POTFAN geometry file is described in this section. The 
current POTGEM version cannot determine all of this data 
because there has not yet been an actual need for it. It 
should be stressed, however, that other POTFAN programs 
assume that geometry files they read in may contain any of 
the data described herein even though POTGEM can not yet put 
it there. 

Data denoted as being default data is only used if not 
changed by other POTFAN programs. The array notation used 
is explained in Appendix C. Variables in the first record 
marked with an asterisk must be entered with a FINISH 
command. All other variables are determined automatically 
or as the result of other commands. Variables or records 
marked with a *■ are ones that would probably not be of 
interest unless POTGEM is used with other POTFAN modules 
(i.e., persons using POTGEM to obtain corner points for 
their own programs may ignore variables and records marked 
with a ♦) . 

In addition to defining the standard geometry file 
format, explanations of how some of the data is determined 
by POTGEM are included. 
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first Hecor d 


(Note--At this point the reader should be familiar with 

Appendices 6 and C.) 

NCTIME — number of words in (CTIME) . 1<NCTIME<5. 

(CTIHE (NCTIHE) ) — Time stamp to identify the approximate 
time of creation of the file. This should be 
printed out in an A4 format whenever the tile is 
read in by a subsequent program. 

NTITL—Nufflber of words in (TITLE). KNTITL<20. 

(TITLE (NTITL) ) — Alphanumeric titling informatxon input 
with the TITLE command. Hhen required, this 
information is to be written out under a format 
such as (1X,20A4). 

NHECS — Number of records on the file includxng the 
first. KNBECS<20. 

(IFOBK (NRECS) ) — An integer describing the format of 
each record. See Appendix B. NRECS and (IFORN) 
are on all POTFAN files so that the EDITPF program 
can be used to list them. In POIGEM the array 
(IFORH) is computed with a FINISH command. 

NID — Number of file identification numbers. 1<NID<10. 

(ID (HID) ) ““Identification number array. In POTGEM 
NID=1 and ID(1) is the same as ID entered with the 
STORE command. Note, however, that if the EDITPF 
or EDITGM program (Section 8.3) is used to modify 
a geometry file created by POTGEM, then NID on the 
modified file will be greater than 1. 

NLOG — Number of words in (LOG) . 1<NLOG<50. 

LOG (1) =flCFLG--This logical variable is .TRUE. iff. 
there are boundary condition flags stbred on a 
subsequent record of the file. For geometry files 
created by POTGEM, BCFLAG will be .TRUE. if a 
BCFLAG command was given or if the null rows of 
panels between geometrical segments have not been 
eliminated by NRI1 and/or NRI2 commands. 

♦ LOG (2) ““Not used any more- 

LOG(3)=OVW — This variable is .TROE. iff, there are unit 
wake vectors in the direction of shed wake lines 
stored on a subsequent record of the file. For 
geometry files created by POTGEM, UVW will be 
•TRUE, if a UVH command was given. 
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♦♦LOG (4) =DANDS--This is .TRUE. iff. both doublet and 
source singularities are to be placed on the body. 

L0G(5)=DSP — This is .TRUE, iff. there are doublet 
singularity flags on a subsequent record. For 
geoaetry files created by POTGEK, DSF uill be 
.TRUE, if a OSFLAG coDiand was given or if the 
null rows of panels between geoaetrical segments 
have not been eliainated by NRI1 and/or NRI2 
comnands. 

L0G(6)=SSF — This is .TRUE. iff. there are source 
singularity flags on a subsequent record. 
Currently POTGEH cannot determine these flags, so 
SSF is always .FALSE, on files created by POTGEH. 

LOG(7)=NTOP — This is -TRUE. iff. there are top surface 
boundary condition vectors stored on a subsequent 
record. If NTOP is .FALSE., then PCTFAN programs 
use the unit normals instead. Currently POTGEH 
cannot determine these vectors, so NTOP is always 
.FALSE, on geometry files created by POTGEH. 

LOG(8)=NBOI — This is -TRUE- iff. there are bottom 
surface boundary condition vectors stored on a 
subsequent record- If NBOT is -FALSE-, then 
POTFAN programs use the negatives of the unit 
normals instead. Currently POTGEH cannot 

determine these vectors, so NBOT is always .FALSE, 
on files created by POTGEH- 

LOG (9) =OTOPL--This is .TRUE. iff- the top surface 
outflow along the top surface boundary condition 
vectors at the control points is stored on a 
subsequent record. If OTOPL is .FALSE., then 
POTFAN programs assume zero outflow. Currently 
POTGEH cannot determine the top surface outflow, 
so OTOPL is always .FALSE. on geometry files 
created by POTGEH. 

LOG(10)=OBOTL--This is -TRUE. iff. the bottom surface 
outflow along the bottom surface boundary 
condition vectors at the control points is stored 
on a subsequent record- If OBOTL is .FALSE., then 
POTFAN programs assume zero outflow- Currently 
POTGEH cannot determine the bottom surface 
outflow, so OBCTL is always -FALSE, on geometry 
files created by POTGEM- 

♦♦LOG (1 1) --Default value for the variable NOSI in the 
VVIM (vortex velocity influence matrix calculator) 
program. Haking the variable -TRUE- if there are 
no bound vortex legs in the N1 direction (e-g.. 
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axisyametric body in axisyanetric flow) will save 
significantly on CPU tiie in VVia. 

♦♦LOG (12) ^-Default value for the variable N0S2 in VVIM. 
Baking this variable .TRUE, if there are no bound 
vortex legs in the N2 direction (e,g., planar 
wings modelled by horseshoe vortices) will save 
significantly on CPU time in VVIM. 

♦♦LOG (13) — Default value for the varicible SLINE1 in 

VVIB. If each of the N2 rows of N1 direction 
vortex legs (i-e. , each group of contiguous N1 
direction panel edges) contains legs that are of 
the same length and parallel to these in its own 
row, then SLINE1 may be set to .TRUE, to save some 
CPU time in VVIB. This situation usually exists 
only for certain, simple, planar components, 

♦♦LOG (14) — Default value for the variable SLIRE2 in 

VVIB. This is the same as SLINE1 except it 

applies in the N2 direction. 

♦ ♦LOG (15) — Default for the variable SUB1 in VVIB. This 

should be .TRUE, if the influence of singularities 
should be summed in the N1 direction (e. g., if the 
component is a shed wake with N1 being the 
streamwise direction) . 

♦ ♦LOG (16) — -Default for the variable SUM2 in VVIB. This 

is the same as SUB1 except it applies in the N2 
direction. 

♦♦LOG (17) =DBLT — This should be .TRUE, if doublet type 
signular ities will be placed on the component. 

♦ ♦LOG ( 18) =SOURCE~This should be .TRUE, if source type 

singularities will be placed on the component. 
Note--LOG(4)=LOG(17) .ARD.L0G(18) . 

♦L0G(19)=RS1 — This iS .TRUE. iff. the force sensing 
locations of the N1 direction vertex segments 
(panel edges of constant doublet singularities) 
are on a subsequent record of the file. If these 
are not on the file, then other POTFAN programs 
use the midpoints determined by averaging corner 
points. These force sensing locations are only 
used if vortices model the component and the 
Kutta-Joukowskii method is used to calculate 
forces, in POTGEB, RSI is established by the PANL 
command. 

♦LOG (20) =RS2--This is the same as BS1 except it applies 
to the K2 direction. 
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NINT — Number of ¥ords in (INT). 1<NINT<50. 

♦INT (1) =NNBC--Nuaber of null boundary condition flags 
(i.e. , number of words in (BCPIAG) that are egual 
to 1) . 


INT(2)=N1 — Total 
direction. 

number 

of corner points 

in 

the 

N 1 

INT (3) =N2 — Total 
direction. 

number 

of corner points 

in 

the 

N2 


INT (4)=N1BC=N 1-1— Number of boundary condition points 
in the N1 direction. This is also the number of 
panels in the N1 direction. 

INT (5) =n 2BC=N2-1— Number of boundary condition points 
or panels in the N2 direction. NIBC^NIBC is thus 
the total number of panels. 

♦INT (6) =NNDS--Number of null doublet singularity flags 
(i.e., number of words in (DSFLAG) that are egual 
to 1) . 

♦ INT (7) =NNS3— Number of null source singularity flags. 

Currently POTGEH cannot compute these flags, so 
NNSS is always 0 on geometry riles created by 
POTGEM. 

INT (8) — Not used. 

INT (9) — Not used. 

♦♦INT(10)=PF1SL— This is the default flow field 

indicator for the 11=1 edge of the component as is 
reguired when vortices model the component. In 
the case of a single segment this is the SL edge. 
A value of 1 is reguired when there will be a 
symmetric image of the component attached at its 
11=1 edge. A symmetric image is one which 

together with the component creates a symmetrical 
flow field about the plane situated symmetrically 
between them. A value of 2 is used for an 
antisymmetrical image. A value of 3 is used when 
the 11=1 edge is physically coincident with the 
I1=N1 edge. Any other value indicates no special 
edge condition. 

♦ ♦INT ( 1 1) =FFISU — This is the same as FFISL ejccept that 

it applies to the I1=N1 edge, 

♦♦INT (12) =FPIVL — This is the same as FFISL except it 
applies at the 12=1 edge. 
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♦♦INT (13) =FFIVU — This is the sane as FFISL except it 
applies at the I2=N2 edge. 

NFLT— Number of words in (FLT) - 1<NFLT<50. 

i 

♦♦FLT (1) =AREF — Component reference area for normalizing 
forces and moments. In POTGEH this is computed 
with a FINISH command if it is not input. See the 
GEOM listing for details. 

♦♦FLT (2) =XLEN 1 — Reference length for normalizing 

momemts about the X axis. In POTGEK this is 
computed with a FINISH command if it is not input. 
See the GEOM listing for details, 

♦ ♦FLT (3) =XLEN2 — Reference length for normalizxng 

moments about the Y axis. In PCTjEH, this is 
computed with a FINISH command if it is not input. 
See the GEOM listing for details. 

♦♦FLT (4) =XLEH3 — Reference length for normalizing 

moments about the Z axis. In PCTGEM, this is 
computed with a FINISH command if it is not input. 
See the GEOM listing for details, 

♦ ♦FLT (5) =DDVax-“-Def ault for the X component of the unit 

waXe vectors. This is used if the wake vectors 
are required and LOG (3) =UVW=. FALSE. 

♦♦FLT (6) =DUVW Y — Same as DUVWX except that it applies in 
the Y direction. 

♦♦FLT (7) =DUVHZ — Same as DUVWX except that it applies in 
the Z direction, 

♦ ♦FLT (8) —Default value for the variable EDVOHT used in 

VVIM. This is the distance beyond which the 
divortlet approximation for a vortex segment’s 
influence will be used. In POTGEM this is 
computed if it is not input. See the GEOM listing 
for details. 

♦ ♦FLT (9) -“Default value for the variable HCUT used in 

VVIM. This is the perpendicular distance within 
which the influence of semiinfinite wake lines is 
set to zero. In POTGEH this is computed if it is 
not input. See the GEOM listing for details. 


♦♦FLT (10) — X 
moments 

component of the 
should be computed. 

position 

about 

which 

♦ ♦FLT(11)— Y 
moments 

component of the 
should be computed. 

position 

about 

which 
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♦*FLT(12)--Z component of the position about which 
moments should be computed. 

S eco nd R ecor d 

(Note — Appendix B defines Jl, J2, J3, and NW.) 

J 1 = N1^-Number of comei: points in the N1 direction. 
J2=N2 — Number of corner points in the N2 direction. 

J3=5 

NW=N1*N2*5 


(XCP(N1,N2) ) 
corner 

— X components 
points. 

of 

position 

vectors 

to 

(YCP (N 1rN2) ) 
corner 

--Y components 

points. 

of 

position 

vectors 

to 

(ZCP(N1,N2) ) 
corner 

— Z components 
points. 

of 

position 

vectors 

to 


♦ (SCP (N1,N2) ) — Values of s at corner points. 

♦ (VCP (N1,N2) ) — Values of V at corner points. 

♦ Nex t Record (if LOG (3) =- 3RUE. ) 

(Note— This record is not present unless 
LOG (3) =-TSUE- .) 

J1 = N1 
J2=N2 
J3=3 ‘ 

HM=N1*N2*3 

(U VBX (N1 ,N2) ) — X components of unit wake vectors along 
wake elements originating from each corner point. 

(OViX (N 1 ,N2) ) — y components of unit wake vectors along 
wake elements originating froo each corner point. 

(OVHZ (N1 f N2) ) —Z components of unit wake vectors along 
wake elements originating from each corner point. 

In POTGEM (DVHX) , (UVHY) , and (UVWZ) are 
established by the DVU command. 
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If this record is not present, then other POTFAN 
prograas assuaethe following: 


UVHX (11,12) 
UVHY (11,12) 
OVBZ (11,12) 


= FLT(5) 

= FLT(6) 

= PLT(7). 


★ N ext R ecor d (if LCG (1) =.TBUE. ) 


J1=N1BC — Number of panels in the N1 direction. 
J2=M2BC — Number of panels in the N2 direction. 


J3=1 

NW=N1BC*N2BC 

(BCFLAG (N1 BC, N2BC) ) — integer boundary condition flag 
denoting the type of boundary condition point. 
Values of 0 indicate a regular boiindiary conditxon 
point. Values of 1 indicate a completely null 
boundary condition point. For closed surfaces to 
be modelled by doublet panels at least one panel 
must have a null boundary condition flag. Any 
other value indicates that the boundary condition 
influence of the panel is to be ignored, but the 
velocity on the panel is to be computed. The 
latter would be the case if constraint functions 
with deleted boundary conditions were to be used. 
If this record is not present, then other POTFAN 
programs assume BCFLAG (II , 12) =0. IN POTGEN this 
array is determined mainly as the result of a 
BCFLAG command. If, however, the null rows of 
panels between geometrical segments are left in by 
not using the NRI1 and/or NE12 commands, then the 
elements of (BCFLAG) corresponding to these panels 
are set to 1, The latter is accomplished with a 
FINISH command and by looking for elements in 
(XBP) (X coordinates of boundary condition control 
points) that are less than or egual to -1.E30. 


♦ Next R ecord (If LOG (5) =. TRUE.) 
JI^NIBC ; 

J2-N2BC 
a 3=1 

N«=N1BC*N2BC 
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(DSFiAG (N1BC , N2BC) ) — Integer doublet singularity ilags. 
These ace used by the VVIH prograu to indicate 
what type of vortex distribution each panel has. 
These flags could also be used by a module that 

replaces VVIM. In VVIM a value of 1 implies no 

singularity; a value of 2 implies a closed 
quadrilateral vortex; and a value of 30 implies a 
horseshoe vortex of type 2 (i.e., a type used for 
wings with 12=1 being the leading edge). There 
are a large nucber of ether types available. The 
various types are shown in figure 6.1-1. Together 
they allow bodies to be modelled with practically 
any combination of vortex segments on panel edges 
in the N1 direction, vortex segments in the N2 
direction, and semi-infinite vortices shed from 
any corner points. If this array is not present, 

VVIM assumes DSfLAG (II ,12) =30. In POTGEM these 

flags are determined in the same ^:ays that 
(BCFLAG) is determined, except that the DSFLAG 
command is used instead of the BCFLAG command. 

♦Kext Re cord (If LOG (6) =. TRUE.) 

J1=N1BC 

J2=N2BC 

J3=1 

NH=N1BC*N2BC 

(SSFLAG (N1BC, N2BC) ) — Integer source singularity flags. 
Currently POTGEM cannot determine these flags. 


♦Mext R ecord 
J1=N1BC 
J2=N2BC 
J3 = 5 

NH = mBCtN2BC^5 

(XBC (NiBC, N2BC) ) — components of position vectors to 
control points. 

(YBC(H1BC,N2BC) ) — Y components of position vectors to 
control points. 

(ZBC (N IBC, M2BC) ) — Z componen ts of position vectors to 
control points. 
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(SBC (NlfiC,N2BC) ) — Values of S at control points. 
(VBC(inBC,N2BC) )— Values of V at Control points. 


♦N ext Hecor d 

d1=NTBC 

J2=N2BG 

J3=4 

NW=N1BC*N2BC*4 

(UNX (N IBC, N2BC) ) — X components of “outward" unit 
normals at control points. 

(UNY (N 1BC, N2BC) ) — Y components of "outward" unit 
normals at control points. 

(UNZ (N 1BC,N2BC) ) — Z components of "outward" unit 
normals at control points. 

The "outward” direction is defined as being on 
that side of the surface on which a vector in the 
N1-cross-N2 direction lies. This may give unit 
vectors that actually point inward. If the 
boun tary condition is one of zero inflow or 
outflow, the direction of the unit normals makes 
no difference. for POTGEW generated geometry 
files, the direction may be reversed by switching 
the SL and SU curves or the VL and VU curves. 

(DA {NlBC,N2BC) ) — Areas of the individual panels. 
POTGEM determines these areas with a FINISH 
command, which in turn calls subroutine ABBAS. 


♦ N ext R eco rd (If LCG (7) =. TRUE. ) 

J1-N1BC 

J2=N2BC 

J3-3 

NW=N1BC*N2BC»3 

(NTOPX (N 1 BC, N2BC) ) — X components of top surface 

boundary condition vectors. 

(NTOPY (NIBC, N2BC) ) — Y components of top surface 

boundary condition vectors. 
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(NTOPZ (M 1 BC,M2BC) — Z components of top surface boundary 
condition vectors. 

The NTOP vectors are not necessarily perpendicular 
to the surface. If this record is absent, 
subsequent POTPAN programs use the unit normals as 
the top surface boundary condition vectors. 

This record is not found on POIGEM generated 
geometry files. 


♦ Next R ecord (if LCG (8) =. TRUE.) 
J1=N1BC 
J2=N2BC 
J3=3 


N«=N1BC*N2BC*3 


(NBOTX (N1BC,N2BC) ) —X components 
boundary condition vectors. 

(NBOTX (N1BC,N2BC) ) — X component: 
boundary condition vectors. 


of bottom surface 


of bottom surface 


(KBOTZ (N1BC, N2BC) ) --Z components of bottom 
boundary condition vectors- 


surf ace 


If this record is absent, subsequent POTFAN 
programs use the negatives of the unit normals as 
the bottom surface boundary condition vectors. 

This record is not found on PCTGEfl generated 
geometry files. 


♦Next R ecor d (if LOG (9) TRUE. ) 

J1=N1BC 

J2=N2BC 

J3=1 

Ni=NlBC*N2BC 

(OTOP (N1BC,N2BC) ) — Desired values of net velocity along 
the top surface boundary condition vectors. If 
this record is not present, then subsequent POTFAN 
programs assume 2 ero net velocity. This record is 
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not found on geometry files created by POTGEM. 

♦Nex t Record (if LOG ( 10) =• IROE. ) 

J 1=N1BC 
J2=N2BC 
J3=1 


NW=N1BC*N2BC 

(OBOT (N 1 BC+N2BC) ) — Desired values of net velocity along 
the bottom surface boundary conditicn vectors. If 
this record is not present, then subsequent POTFAN 
programs assume zero net velocity. This record is 
not found on geometry files created by POTGEM . 


♦Next Re cor d 


This record and the following three records define the 
outside boundaries of the component. This definition 
is necessary because panels need not extend up to the 
boundaries. 

J1 = N1 

J2^10 

J3=1 


NW=N1*10 

(XVLC(N1) — X values at the intersection of the 11=1 
boundary (VL(S) curve in POTGEM) with the corner 
point grid lines extending nominally in the V or 
N2 direction. If the panels extend to the 
boundaries of the component, then 

XVLC (12) =XCP ( 1 , 12} , but panels need not extend to 
the component boundaries (e.g,, the panels at 
leading edges cf thin wings to be modelled with 
voftiees are set back from the edge). 

The significance of most of the arrays on this and 
the next three records should be clear from the 
explanation of (XVLC), therefore, they will not be 
defined. 

(XVUC(NI)) 

(YVLC(NI)) 

(YVOC(NI)) 

(ZVLC (HI) ) 


6-12 



(ZVUC (N1) ) 
(5VLC (N1) ) 
(SVUC (N1) ) 
(VVLC (Nl) ) 
(VVUC (Nl) ) 


♦Next Recor d 
J 1=N1BC 


J2=12 


NH=N1BC»12 

(XVLB(NIBC)) — Same as (XVLC) except that boundary 
condition point grid lines are involved. 

(XVUB(NIBC)) 

(YVLB(NIBC) ) 

(YVUB(NIBC)) 

(ZVLB(NIBC)) 

(ZVUB(N IBC) ) 

(SVLB(NIBC) ) 

(SVUB(NIBC)) 

(VVLB(NIBC)) 

(VVUfl(NlBC)) 

(CORD1 (NIBC) ) — Reference chord lengths in the N2 
direction. These are used in PCTFOB, which is the 
program that computes span loads. In POTGEM this 
array is computed after a FlhlSH command. 

(SPANi (N1BC) ) — Reference widths of rows of 

singularities extending in the N2 direction. 
These widths are some measure of lengths in the Nl 
direction and are used in determining sectional 
aerodynamic properties such as section lift 
coefficient- In POTGEM this array is computed 
after a FINISH command. 


♦Next R ecor d 
J1=N2 


J3=1 

NB=N2*10 
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{XSLC (N2) ) 
(XSUC (N2) ) 
(YSLC (N2) ) 
(YSrJC (N2) ) 
(ZSI.C (N2) ) 
(ZSUC (N2) ) 
(SSLC (N2) ) 
(SSUC (N2) ) 
(VSLC (N2) ) 
(VSUC (N2) ) 


♦N ext Record 
J 1=N2BC 
J2=12 


J3=1 

NW=N2BC*12 


(XSLB(N2BC)) 

(XSUB(N2BC) ) 

(YSLB (N2BC) ) 

(YSUB(N2BC)) 

(ZSLB(N2BC)) 

(ZSUB (N2BC) ) 

(SSLB(H2BC)) 
iSSUB(N2BC) ) 

(VSLC(N2BC)) 

(VS0B(N2BC)) 

(CORD1 {K2BC) ) 

(SPAN2 (N2BC) ) 

♦ Next Rec ord (if LOG { 19) =.TRUE. ) 

J1=N1BC 

J2=N2 

J3=3 

NH=N1BC*N2*3 

(XS 1 (N1BC,N2) ) — X components of force sensing locations 
on the N1 direction vortex segments. 

(YSl (N 1BC,N2) ) --Y components of force sensing locations 
on the N1 direction vortex segments. 

(ZS 1 (N 1BC, N2) ).— Z components of force sensing locations 
on the N1 direction vortex segments. 
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These arrays are used only it vortices are used to 
model the component and then only if the Kutta 
Joukowskii theorem is used to determine the forces. If 
this record is absent, subsequent POTFAN programs 
requiring the N1 force sensing locations assume the 
midpoints of the N1 segments. In POTGEM these arrays 
ace computed if the variable RSI is entered as .TRUE, 
on a PAML command. 

±N§xt_R^ord (if LOG (20) =. TRUE.) 

J1 = N1 

J2=N2BC 

J3 = 3 

NW=N1*N2BC*3 

(XS2 (N1 ,R2BC) ) — X components of force sensing 
locations on the N2 direction vertex segments. 

(Y S2 (N 1 ,N2BC) ) --Y components of force sensing locations 
on the N2 direction vortex segments. 

(ZS2 (N 1 , N2BC) ) — Z components of force sensing locations 
on the N2 direction vortex segments. 

These arrays serve the same function for the N2 vortex 
segments that (XS1) , (YSl), and (ZS1) do for the N1 
vortex segments. In POTGEM these arrays are calculated 
if the variable RS2 is entered as .TRUE, on a 
PARL=COMMAND 


6.2 SUMMARY OF GEOMETRY FILE DATA 

This section summarizes the geometry data by presenting 
in pseudo-FORTRAN form the statements that could be used to 
create a geometry file. They are as follows: 

CALL OPENW(NTG,1,ID(NED),1) 

WHITE (NTG) RCTIME, (CTIME) ,NTITL, (TITL) , 

#NHBCS, (IFORM) ,RIE, (ID) , 

#NLOG, (LOG) ,NINT, (INT) ,NFLT, (FLT) 

R1=INT(2) 

N2=INT(3) 

N1BC=INT(4) 

N2BC=INT (5) 

WRITE (NTG) Nl,N2,5,Hl*N2+5, (XCP) , (YCP) , (ZCP) , (SCP) , (VCP) 
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IF (LOG (3)) WHITE (NTG) N 1 , N3 ,2, K 1* N2^3 , (OV WX) , (U VW X ) , ( UVWZ) 
IF (LOG(1)) WRITE (NTG) N 1 BC, N2BC, 1 , Hi BC* N2BC, (BCFL AG) 

IF (LOG(5)) WRITE (HTG) N 1 BC, N2BC , 1 , N 1 BC*N2 BC, (DSFLAG) 

IF (LOG (6)) WRITE (NTG) N 1BC, N2BC^ 1 , N1BC*N2BC, (SSFLAG) 

WRITE (NTG) N1BC,N2BC,5,N1BC*IJ1BC*5, 

# (XBC) , (YBC) , (ZBC) , (SBC) , (VBC) 

WRITE (NTG) N1BC,N2BC,4,N1BC»H2BC*4 

# (UNX) , (UNY) , (UMZ) , (DA) 

IF (LOG(7)) WRITE (NTG) N1BC,N2BC,3, NlBC*N2BC+3, 

# (NTOPX) , (NTOPY) , (NTOPZ) 

IF (LOG (8) ) WRITE (NTG) N 1 BC, N 2BC, 3 , N 1BC* N2BC*3, 

# (NBOTX) , (NBOTY) , (NBOTZ) 

IF (LOG (9)) WRITE (NTG) N1BC,N2BC, 1, N1BC*N2BC, (OTOP) 

IF (LOG (10)) WHITE (NTG) N1BC,N2BC,1 ,N1BC*H2BC, (OBOT) 

WRITE (NTG) N1, 10, 1,N1*10, (XVIC) , (XV DC) , (Y VLC) , (YVUC) , 

# (ZVLC) , (ZVUC) , (SVLC) , (SVUG) , (VVLC) , (VVUC) 

WRITE (NTG) N1BC, 12, 1,N1BC*12, (XVLB) , (XVUB) , (YVLB) , (YVUB) , 

# (ZVLB) , (ZVUB) , (SVLB) , (SVOB) , (VVLB) , (VVUB), (CORD2) , (SPANI) 

WRITE (NTG) N2, 10, 1,N2*10, (XSLC) , (XSUC) , (YSLC) , (YSUC) , 

# (ZSLC) , (ZSUC) , (SSLC) , (SSUC) , (VSLC) , (VSDC) 

WHITE (NTG) N2BC, 12, 1,N2BC*12, (XSLB) , (XSUB) , (YSLB) , (YSUB) , 

# (ZSLB) , (ZSUB) , (SSLB) , (SSUB) , (VSLB) , (VSUB), (CORDI), (SPAN2) 

IF (LOG (19)) WRITE (NTG) N1BC,N2,3,N1BG*N2*3, 

I (XS1) , (YS1) , (ESI) 

IF (LOG (20) ) WRITE (NTG) N1,N2BC,3,N1*N2BC*3, 

# (XS2) , (YS2) , (ZS2) 

In POTGEM the geometry file is written by subroutine 
STORGJl. 
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7 SAMPLE CASES 


This section presents a number of sample cases. These 
cases were not only devised to test the program, but to 
serve as the basis of a tutorial guide for assisting users 
in learning how to use the program. This tutorial proceeds 
from simple cases to more complex ones with each new facet 
of the program being explained the first time it is 
introduced. Potential users are therefore advised to study 
every sample case. 


7. 1 TES T CASE NO. 1 - T H IN, SYHHETBICAL. SWEPT. F LA T liilNG 

The so-called Warren 12 planfora shown in Figure 7.1-1 
has been chosen as the first test case. Note the SL(V), 
SU (V) , VL (S) , and VU{S) curves. This wing consists of a 
single segment, which is the right hand half of the wing, 
and only half the wing needs to be considered due to the 
symmetry- The spanwise direction has been chosen as the S 
or N1 direction and the chordwise direction as the V or N2 
direction. These could have been reversed, but for more 
complicated wings this is generally the best way to choose 
the directions. This wing can be handled most simply by 
using the WING command or else by inputting the SL, SU , VI, 
and VU curves explicitly. The former method is the subject 
of this section, while the latter is discussed in Section 
7.2. 

The input deck illustrating the use of the WING command 
is shown in Figure 7-1-2. Comments concerning some cards 
are given below: 

Card 7 — This card specifies that there should be 10 panels 
in the S direction (spanwise in this case) and 4 in the V 
direction (chordwise) . 

C ard 9 — Uniform spacing of the vortices in the S direction 
is requested. 

C ard 1 1 — IOPT=2 specifies the typical (1/4 and 3/4) spacing 
of vortices and control points in the V direction. 

Ca rd 13 — BS1 is -TRUE. because it is intended that loads 
should be determined using the Kutta-Joukowskii law. Ther 
will be no load-carrying vortex segments in the N 
direction. 
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Ca rds 14-1 7 — These cards cause all panels to have a type 30 
doublet singularity on then. This type of singularity is a 
horseshoe vottex with bound portion in the N1 direction. 
The -1 values on cards 15 and 16 signal the program to use 
the number of panels in the N1 and N2 directions, 
respectively, as the upper limits. 

C ard 19 — LOG (12) is .TEUE. because there are no bound vortex 
legs in the N2 direction. 

LOG (13) is .TRUE. because the vertex legs along each panel 
line in the N1 direction are all parallel and of the same 
length. 

IliT(IO) is 1 because the SL edge (i.e., the centerline) is 
in a symmetry plane. 

FLT(5), FLT (6) , and FLT (7) are the components of the unit 
wake vectors. In this case the wake vectors are all 

identical (if the angles vary from corner point to cornet 
point the UVW command should be used) and parallel to the X 
AXIS- 

The printed output for this case is shown in Figure 
7.1-3. Figure 7.1-4 shows the arrangement of panels, wing 
outline (dashed), and wake vectors. This figure was 

generated using the PLOTGM program (Section 8.4). The input 
to PLOTGM is shown in Figure 7.1-5. 


7.2 TEST CASE MO. 2 

The results from this test case ate identical to case 
no» 1. The only difference is that this case was run using 
a different input method. With this method, which is more 
complicated, but more general, it is necessary to explicitly 
consider the component axis, the cross sectional planes, and 
the VL, VU, SL, and SU curves. The component's axis is 
taken to be the Y axis. The Y» axis is the X axis and the 
Z* axis is the Z axis. Thus the Y-Z plane must be rotated 
by -90 degrees about the Z axis in order to become parallel 
to the Y'-Z* plane. The SL, SU, VL, and VU curves are shown 
in Figure 7,1-1. The input is shewn in Figure 7,2-1. 
Specific comments on the input deck follow: 

Card 6 — COPT=0 causes the corresponding function to be zero. 
Thus XAX1S(S)=0. 

Card 8 — COPT=2 implies YAXIS(S)=S. 

Ca rd 10 — Causes PHI(S) to equal -90. 

C ard 1 2 — Causes EX (S) to equal 0. 
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C ard 14 — Causes EY(S) to equal 1. Note that NTAB =1 does 
not have to be entered because it was given on Card 10. 

C ayd_ 16 — This is the only cross section. Since there is 
only one cross section, SCS does not have to be entered. 
COPT- 0 then implies that V2(S,V)=0 or Z(Y,X)=0. 

Cards 23-30 — These define the Si, SU, VL, and YU curves 
shown in Pigure 7.1-1. 

The output from this case is shown in Figure 7.2-2. 

This sane wing could have been done by a sonewhat 
different procedure. Nanely the default axis and cross 
sections could have been used initially and, after 
panelling, the entire conponent could have been rotated by 
90 degrees and about the Z axis. If this were to be done, 
cards 5 through 14 would be deleted and a ROSS command and 
appropriate data would be inserted after card 32 or 37. 

7.3 T EST CASE HQ. 3 NA S A ANE S 12.' 1 92-by-24.384 ttETER WIN D 
TU NNEL 

A cross section view of the right half of the 
12. 192-by-24. 384 (40*x80») test section is shown in Figure 
7.3-1. It can be described analytically as follows; 

Section A-B; Z*(7*) = -6-096 for -6.096<Y»<0 

Section B-C: r (©) = 6.096 for -90 degrees < 6 <90 degrees 

Section C-D: Z»(y*) = 6.096 for 0 >Y* > 06.096. 

The most accurate method of describing this geometry to 

the program is to take advantage of the multiple segment 
capability with the segments being the sections A-B, B-C, 
and C-D. The V Variable is y* in sections A-B and C-D and 
is e in section B-C. The S variable is chosen to be X. The 
equations for the axis are the following: 

XAXIS (S) - S 
YAXIS (S) - 6. 096 
ZAXIS(S) = 0. 

In the first segment (A-fi) the equations for V2, VL(S}, and 
VD(S) are: 

V2 = -6.096 
Y1 = -6.096 
YU = 0 

In the second sega.ent (B-C) the equations for V2, VL(S|« and 
VU(S) are: 
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V2 - 6.096 

VL = -90 
VU = 90 

In the third segnent (C-B) the equations for V2, VL(S), and 
VU(S) are: 


V2 = 6.096 

VL = 0 

VU = -6.096 


Note that VU (S) 


may be less than VL(S). 


If the longitudinal extent of the portion of the tunnel 
to be modelled is from X = -10 to X = *20, then SL(V) and 
SU(V) are: 


SL = -10 
SO = 20 

The input which implements the description just given 
is shown in Figure 7.3-2. Comments on certain of the input 
cards are given below: 

Cards 4-5 — Define XAXIS(S) as 6.0 96. 

C ard 6 --Defines V2 as y*. 

Ca rd 1 0 — NSEGVT is the total number of segments in the V 
direction. NBPV - 3,10,3 causes 3 panels in the first 

segment (A-B) , 10 panels in the second segment (B-C) , and 3 
panels in the third segment (C-D) . 

Ca rd 23 — Panels the first segment. 

Ca rd 25 --Defines V2 as 6. 

Ca rds 2 8-29 — Name the next segment to be panelled. 

Cards 36-37 — Panel the second segment. Note that the VLBC 
command was not required for this segment. 

Cards 47-48 — Panel the third segment. Note that the VLBC 
command was not required for this segment and also it was 
not necessary to give ah SBIl command for IC=11 because V2 
is numerically the same as in the previous segment (even 
though its basic definition has been changed by the CARY 
command on card 38). 

Cards 5 3-56 — Set all doublet singularity flags to type 2. 

Cards 57 -59 — Set the doublet singularity flags along the 
trailing edge to type 4. See Figure 6.1-1. 
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The output from the program is shown ia Figure 7.3-3- 
Figure 7.3-4 shows various views of the resultant panel 
distributions, control points and shed wakes- Figure 7,3-5 
shows the PLOTGM input that generated Figure 7-3-4. 


7 . 4 TEST CASE NO. 4 - THIN ^ S WEP T, O MCAHBERED , UNTWISTED 

iiiSISiillEiiiMiAL " ’ 

Top and rear views of the subject wing are shown in 
Figure 7.4-1- 

The first step in treating this wing is to rotate it by 
90 degrees around the Z axis. The resultant plan view is 
shown in Figure 7.4-2. Note the Si, SO, VL and VU curves 
and that S = -X and V = y* = Y. Also V2 = z* = Z. Two 
cross sections which are sufficient to complete the geometry 
description are: 

V2-z»= 0atS=0 

and 

V2 = z* = .5 at S = 2 

After the wing is input to the program via the above 
description, it is rotated by 90 degrees about the -Z axis 
so that it conforms to Figure 7.4-1 and then it is rotated 
by 45 degx:ees about the Y axis to put it at 4 5 degrees angle 
of attack.. This final rotation is not necessary to obtain a 
solution for the wing at 45 degrees angle of attack, but has 
been done for testing purposes. 

A schematic of the vortex singularity model used for 
this wing together with the doublet singularity flags is 
shown in Figure 7,4-3. 

The input deck for this case is shown in Figure 7.4-4. 
Comments on some cards are given below: 

C ard 6 — COPT=3 defines XAXIS (S) as being equal to -S. 

Cards 7-10 — These cards define the two cross sections. 

Ca rds 2 8-29 — These cards rotate the wing so that the YAXIS 
is the spanwise direction. 

Cards 3 2 - 42 - — These cards define the doublet singularities 
shown in Figure 7.4-3. 

Card 44 — LOG(13) is -TRUE, because the N1 vortex segments 
for any given value of N2 are all of the sane length and 
parallel. 

The output from the program is shown in Figure 7.4-5. 
Figure 7.4-6 shows top and side views of the wing outline 
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(dashed) , panels, and shed wakes 


7.5 T EST CASE NO. 5 - S PHE RE WITH S THE CIH C UMf ERENTIAL 

VAR IABL E 

Figure 7.5-1 shows the choice of the S and V variables 
and the y*-z* plane, which varies with S. Note that as S 
varies, the origin of the y*-z» plane remains fixed. 
Therefore, it is necessary that 

lAXIS(S) = YAXIS(S) = ZAXIS(S) = 0. 

The variable PHI in this case is identical to S and the 
-Y axis is the rotation axis Q’l the cross section. Hence, 

PHI(S) = S 
EX(S) = 0 
EY (S) = - 1 
EZ(S) - 0. 

Also the variable V2 in this case is just the radius of the 
sphere and is taken egual to 1. 

Figure 7.5-2 shows the input deck for this case. A 
total of 5 panels are laid out in the S direction with S 
varying from 1 to 90 degrees (i.e., only 1/4 of the sphere 
is panelled) . The panel corners occur at egual increments 
in S. In the V direction 10 panels are laid out at egual 
increments with the first cornet points being moved back 
from the nose (V = 0), which is necessary for the sphere to 
be modelled with doublet singularities. Type 19 and type 27 
doublet singularities are selected due to the axisymmetry 
(see Figure 6.1-1). Finally the component is rotated so 
that the nose is on the negative X axis. 

Figure 7.5-3 shows the output from POTGEH for this 

case. 

7 . 6 TEST CAS E NO. 6 - T WO DIMENSIONA L A IRF OIL 

The purpose of this test case is to show how cross 

sections of typical thick wings can be handled. The cross 

section is decomposed into two segments. The first segment 
is the lower surface starting at the trailing edge and 
ending at the nose. The second segment is the upper surface 
beginning at the nose and ending at the trailing edge. For 
loth cases the V- variable is X. 

In both the first and second segments the input data is 
proportional to S,QRT(V) near V = 0. The interpolation 

routines provided in the program are not able to fit curves 

to such data accurately. To overcome this difficulty one of 
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the transformation capabilities of the prograi is invoked so 
that the independent variable for interpolation is (V) . 
In terns of this variable the input data is nov linear near 
V = 0, and, therefore, interpolation is much more accurate. 
See Section 3, 1-4 for more details- 

The program deals with three-dimensional bodies, 
therefore the airfoil section is considered to be of unit 
thickness and centered at the origin. Thus the S variable 
is Y and SL(V) = -.5 and SU(V) = -5. Also there is one 
panel in the S direction. 

The section is input as if it were parallel to the Y-Z 
plane and after the section is panelled, it is rotated so 
that it is parallel to the X-Z plane. This is the typical 
procedure that is followed for all wings when the usual 
aerodynamicist *s coordinate system is used. 

The input deck for this case is shown in Figure 7.6-1 
and the program output is shown in Figure 7.6-2. Specific 
comments on the input deck arc given below: 

C ard 6 — COPT=6 specifies that the controlled deviation 
interpolation method be used for interpolation within the 
cross section. IOPT1=1 specifies that the independent 
variable for this interpolation will be SQBT(V) instead of 

V- 

Card 29 — IOPT=6 gives dense panels near the nose- 

Cards 7 3-80 — When modelling closed bodies with doublet 
panels, a hole must be left somewhere in the body. These 
cards make the first panel null and thus the hole will be at 
the trailing edge on the bottom surface. These cards also 
define the remaining panels as type 32 (see Fig. 6.1-1). 

C ard s 81-88 — These cards create a null boundary condition 
point” at ”the null panel so that the resulting influence 
matrices will be square. 

Ca rds 89 -9 5 — In order to make the flew field two 
dimensional, the wake vectors on the II = 1 edge must be 
directed in the +Y direction while the vectors on the II = 2 
edge must be directed in the -Y direction. 

Figure 7.6-3 shows the resulting corner point model of 
this airfoil. 


7 . 7 TEST C ASE HQ. 7 - T HIN WIMG WITH TWIST , CAtlBEfl, AN D 
DIHEDiAl 

The purpose of this test case is to show how 
complicated, thin wings can be handled. 
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The planfora of the subject wing prior to the addition 
cf the twist is shown in Figure 7.7-1. At ¥ = 0.0 the uean 
camber line is Z/c =0.0 where c is the wing chord- At Y = 
15.479 the camber line is the parabolic arc given by 

Z/c = X (1-x) 

where x is the linear function of X that is equal to 0.0 at 
the leading edge and 1.0 at the trailing edge. At Y = 25.0 
the camber line is derived from Table 7.7-1 by multiplying 
the Z/c values by a factor of 12.0. Note that the table 
entries extend beyond the leading edge. This is done in 
order to control the interpolation near the leading edge. 
Also it should be noted that between X = .2025 and X = 1. 0 
the camber line is exactly straight. At Y = 35.107 the 

camber line is also derived from Table 7.7-1, only the 
factor to be applied is 15.0. Between the given stations 
the mean camber surface varies linearly with Y along the 
lines X = constant. 

The twist in degrees is equal to -.0097363* Y**2. The 
twist center is the trailing edge of the wing and the twist 
axis is everywhere parallel to the Y axis. The vertical 
displacement of the trailing edge of the wing is given by 
Z = .25Y. 

If the cross sections are initally chosen to be 
parallel to the Y-Z plane, the above description of the 
geometry can be accomodated by the following: 

S = X 

YAXIS(S) = -Xte(Y) (Fig. 7.7-1) 

ZAXIS(S) = .25*3 

PHI(S) = -.0097636*S**2 
EX(S) = 1 
EY(S) = 0 
YPSCAL(S) = c (Y) 

V = 1 - X = (Xte - X)/c(Y) 

VL(S) =1 
VU(S) = 0 

In addition V2 will be the value of Z on the surface prio r 
to adding in the twist or dihedral. 

The wing is modelled with quadrilateral vortices. 
These are placed in the typical manner chordwise and in both 
the root segment and tip segment they are evenly spaced 
spanwise. A waXe is shed from the trailing edge and the 
rearward half of the tip edge and trails back at an angle of 
10 degrees. 

This wing is treated in a manner similar to test case 
no. 4. The primary differences are the addition of camber 
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and twist, the dihedral is handled soaewhat differently, and 
the presence of the break in the trailing edge creates the 
need to assure that panel edges occur also at the break. 

The input deck for this case is shown in Figure 7.7-2. 
Specific comments on some of the cards are given below: 

Ca r d 35 — C0PT=6 and PAHAM(1)=0 cause the controlled 
deviation interpolation method with linear interpolation in 
the end intervals to be used for this cross section. Linear 
interpolation in the last interval is ideal because the rear 
portion of the camber line is exactly straight. Linear 
interpolation near the leading edge would not be very good, 
however, so the table has been extended beyond the leading 
edge. 

Ca rd 42 — The table (Cards 36-39) does not have to be entered 
again even though it was transformed by the previous 
command. This is because the transformations are not done 
cn the intermediate arrays (see subroutine SRFIhl). 

Ca rd 68 — Note that the SL (V) curve does not have to be 
entered for the second segment. This is because it is 
identical to the SU(V) curve from the previous segment. 

Car ds 6 1, 84, 8 7, a nd 9 0 — See Figure 6.6-1 for a schematic 

of these vortex models. 

The output from this case is shown in Figure 7.7-3. 
Various views of the panels, boundary condition points, and 
shed wakes are shown in Figure 7.7-4. 
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8 ttELATtP PR0GHAM 5 

This chapter describes prograas and subrcutines related 
to POTGEH and which are also available. 


8-1 S UBROUTI NE BEA DGM 

This subroutine reads a geometry file in a 
straxgh tf orvard manner. The principal drawbacks are that 
all of the data on the file is read in whether necessary or 
not and the program does not supply the data defaults 
discussed in Section 6.1. 


8.2 SUBRO UT INE HDGHA 

This subroutine also reads in a geometry file. It 
differs from EEADGM on three counts, first cf all, it does 
not read in any data that is not required. Secondly, if 
some data is not available on the file, but is requested, 
BDGNA will compute it according to established conventions 
(see Section 6.1). Finally BDGHA packs all of the requested 
data solidly into a single array and computes the addresses 
of the individual arrays. Thus FDGMA is compatible with 
calling programs using dynamic memory allocation. 


8.3 EDITGM 

EDITGM is a separate program that can be used to 
manipulate POTFAN geometry files in various ways. It 
operates in the same way that all POTFAN programs do. That 
is, it operates using commands as the basic form of logic 
control. The commands in EDITGM allow geometry files to be 
read in, printed out, edited, compared with other geometry 
tiles, and, after editing, to be stored. Also, EDITGM can 
be used to rotate, shift, and scale the component. 

The EDIT command has been found to be very useful for 
changing flags ((BCFLAG) and (DSFLAG) } and parts of the 
introductory record that should have been entered with a 
FINISH command but were not. 

The PRINT command is another useful command. It is 
frequently used to print out geometry file data that, for 
whatever reason, was not printed cut when POTGEM was run- 
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The BOSS command can be used to rotate, shift, and 
scale the component. This command works in the identical 
way that HOSS works in POTGEM. 

The listing of EDITGM contains all of the necessary 
instructions to run the program. 


8.4 PlOTGil 

PLOTGH is an independent and separate program that can 
be used to graphically display a component described by any 
FOTFAN geometry file. This display consists mainly of the 
corner points connected by straight line segments in either 
or both directions. Optionally the boundary condition 
points, unit normals, unit wake vectors, and component 
outline can be displayed. The display can be from any 
viewing angle and to any desired scale. 

PLOTGa is currently set up to run a Zeta plotter, but 
can be easily converted to use a Calcomp plotter. 

PLOTGM works on a command basis in the same way that 
all other POTFAN programs do. These commands and all 

necessary user instructions are contained in the PLOTGM 
listing. PLOTGM was used to prepare figures for this 

report. In some cases the input decks have been shown to 

serve as examples since no separate documentation will be 

available for this program (see chapter 7). 
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APPENDIX A 


A— STANDAapIZED FI LE HANDLI NG PRO CEDD RES FO B POTfAN PRO GEA HS 


Standardized FOBTfiAN procedures and subroutines tor 
opening and closing tiles have been developed to facilitate 
using and coding POTFAN programs and the conversion of these 
codes to different computet systems. 


A. 1 FILE CREATION 


This section describes actions taken before and after 
any POTFAN program attempts to write a POTFAN file. 

Prior to writing any permanent file onto a unit, all 
POTFAN programs call a system dependent subroutine as 
follows : 


CALL OPENW (NT,IFTYP,ID,IR) 


If IR is not zero, then NT and ID are considered subroutine 
inputs- NT is the logical unit number on which the file 
will be written and ID is the file creation identifier, 
which should also be the primary file identification number. 
If IB is zero, then ID is net considered a subroutine input 
and NT is only the default unit number. In this case the 
program reads in ID and NT from a card via 215 format. If 
the value of NT on the card is zero, the subroutine replaces 
NT with the default value. 

If the value of ID determined in either case is then 
still zero and if it is possible on the computer system 
being used, the program will replace ID with the current 
number on the identification number file and also update the 
identification number file. 

In addition to NT, ID, and IB, IFTYP is also input to 
the program. IFTYP defines the type of file being created 
according to the following table: 



IFTYP TYPE OF FIL E 


1 Geometry 

2 Boundary condition 

3 Influence matrix 

4 Velocity matrix 

5 Solutions 

6 Velocity at force sensing location of N1 segments 

7 Velocity at force sensing location of N2 segments 

8 Constraint function transformation matrix 

9 Zeta plot file 

10 Constrained influence matrix 

12 Preset solution 

15 Inverse or decomposition of influence matrix 

16 External velocity 

17 Surface pressures 

18 Surface velocity 


Once ID and NT have been determined, the program opens 
(if possible on the system being used) the file for writing 
using a file name determined from ID and IFTYP. On IBM 
systems, opening a file consists of issuing a DDEF to the 
operating system. On the INFONEI UNIVAC 1108 system, an 
EQUATE command is involved. This feature eliminates the 
need for job control cards to handle files on those systems 
for which FORTRAN programs can open files. 

The program then rewinds the file and writes a message 
indicating which unit has been opened and the value of ID 
and IFTYP. 

After the file has been opened and written upon, it is 
released by calling another system dependent subroutine as 
follows: 


CALL ENDFIL (NT) 


This subroutine writes an end-of-file mark on the unit 
and (if required by the system being used) , releases the 
unit. The subroutine also writes a message indicating that 
unit NT has been closed. 
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A. 2 FILE ACCESSING 


This section describes actions taken before and after 
any POTFAN program attempts to read any POTFAN file. 

Prior to reading any permanent file from a unit all 
POTFAN programs call a system dependent subroutine as 
follows : 


CALL OPENB(NT,IFTYP,IO,IB) 


It IR is not zero, then NT and ID are considered 
subroutine inputs. NT is the logical unit number from which 
the file is read and ID is the file access identifier, which 
should also be the primary file identification number. If 
IR is zero, then ID is not considered a subroutine input and 
NT is only the default unit number. In this case, the 
program reads in ID and NT from a card via 215 format. If 
the value of NT on the card is zero, the subroutine replaces 
NT with the default value. 

In addition to NT, ID, and IB, IFTYP is also input to 
the program. IFTYP defines the type of file being read 
according to the table in the previous section. 

Once ID and NT have been determined, the program 
attempts to open the file using a file name determined from 
ID and IFTYP. The capability to open a file from a FORTRAN 
program depends on the system being used. As explained in 
the previous section, this may involve a DDEF or EQUATE 
command and can eliminate the need for job control cards to 
handle files. 

The program rewinds the file and writes a message 
indicating which unit has been opened and the value of ID 
and IFTYP. 

After control is returned to the calling program and 
the first record of the file has been read, all POTFAN 
programs check to see if the access identifier is equal to 
the actual primary file identification number existing on 
the first record. If not equal, the program writes an 
informational diagnostic message and proceeds. This feature 
is meant to be a helpful filekeeping technique for those 
systems that do not permit automated file control. 

After the file has been read and there is no further 
use for it, it is released by calling another system 
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dependent subroutine as follows: 


can FILEND (NT) 


This subroutine rewinds unit NT and (if required by the 
system being used) releases the unit. 
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APPENDIX B 


E ST ANDARDIZED FOR MAT O F POTF AS FILE S 


A standard fornat has been developed for POTFAN files. 
This format is applicable to all files except scratch files 
and plot files. This standard has been developed for the 
following reasons: 


1. to minimize the effects of changes in one POTFAN 
segment on other POTFAN segments; 

2. to allow a program to be developed (EDITPF) which 
can list and/or edit the contents of any POTFAN file; 
and 

3. to promote consistency among POTFAN programs. 


Briefly, the standardized POTFAN file consist of one or 
more records. The first record is called the introductory 
record and contains miscellaneous data including the primary 
identification number, a title, and real, integer, and 
logical parameters reflecting how the data on the remainxng 
records was calculated and/or how it is to be used. The 
second and subseguent records generally contain the bulh of 
the data and are called data records. The latter records 
contain one or more arrays which are always either integer 
cr floating point numbers (i.e. integer and floating point 
numbers are not mixed on a single record) . A detailed 
description is given below. 


F irst R eco rd (Introductory Record) 


This record is created by an unformatted write 
statement such as the following: 

WRITE (NT) NCTIME, (CTIME (N) , N=1 , NGTIKE) , NTITL, 
i (TITI. (N) ,N=1 ,NTITL) ,NRECS, (IFORH (N) ,N= 1 ,NRECS) , 

INID, (1D (N) ,N=1,NID) , NLOG, (LOG (N) , N= 1 , NLOG) , 
iNINT, (INT (N) ,N=1,NINT) , NFLT, (FLT (N) , N=1 , NFLT) 



The values of NCTIflE, NRECS, KID, NLOG, NINT, and NFLT 
are all at least one and can vary from file tc file even for 
files of the same type (e.g. NINT aay be different on two 
different geonetry files) . An explanation of these 
variables is given below: 


NCTIME Number of words in (CTIME) 


(CTIhE) 


NTIIL 

(TITL) 

NBECS 

(IFORB) 

BID 

(ID) 

NLOG 

(LOG) 

NINT 

(INT) 

BELT 

(FLT) 


Creation time in A4 alphanumeric format. Whether 
or not this array can be filled out depends on the 
availability of a system dependent subroutine to 
compute it. This array is used only as a 
filekeeping aid. It is printed cut whenever a 
file is created or read. 

The number of words in (TITL) . Generally NTITL is 
a multiple of 20. 

Alphanumeric titling information (e.g. "Delta wing 
with flaps"). This array is to.be written under a 
format such as (1x,20A4/). 

The number of records (including the first) 
comprising the file. NBECS is also the number of 
words in (IFORM). 

An integer array indicating the kind of numbers on 
each record. A value of zero implies an integer 
and a value of one implies a floating point 
number. IFORM (1) has no significance. 

The number of words in (ID) 

Identification number array. IC(NID) is the 
primary file identification number. In order to 
keep track of files ID (NID) should be unique for 
each file. This number is printed out whenever 
the file is created or read. 

The number of words in (LOG) 

An array of logical parameters 

Number of words in (INT) 

An array of integer parameters 

Number of words in (FLT) 

An array, of floating point parameters. If the 
remaining data on the file is dependent on Mach 
number, then FLT(1) is the Mach number. 


B-2 



S econd and S ub s equen t R ec ord s (Data Records) 


The renaining records o£ POTFAM files contain one or 
■ore arrays. If the data record contains more than one 
array, then all arrays on the record must be of the same 
type (i.e. either integers or real numbers, but not both) 
and all arrays must have the same number of words. The 
records also contain array dimensions (J1, J2, and J3) and 
the total number of words in all arrays on the record (NU) . 
Following are some examples of code used to create data 
records: 


NW = J1*J2*J3 

WaiTE(NT) J1,J2,J3,NM, ({(A(I,J,K) ,I=1,J1) ,J=1,J2) ,K=1,J3) 


J3 = 2 

NH - J1»J2«J3 

WRITE (NT) J1,J2,J3,NW, ( (A (I,J) ,1=1 ,J1) ,J=1 ,J2) , 
#((B(1,J) ,I=1,J1),J=1,J2) 


J2 = V 
J3 = 1 
NU = J1 

WaiTE(NT) J1,J2,J3,NK, (A (I) ,1= 1, N« ) 


J2 = 3 
J3 =1 
NH = 3*J1 

WRITE (NT) ai,J2,J3,NH, (A(I) ,I=1,J1), (B (I) ,I=1,J1) , 
#(C(I) ,I=1,J1) 


Note that in the above examples all dimensions with 
■ultiple arrays were written with the leftmost indices 
varying most rapidly. This practice is always followed 
unless it is strictly necessary to do otherwise. 

No matter how a data record was created, it can be read 
in by either of the following: 


BEAD (NT) J1,J2, J3,NN, (A (I) ,1=1 ,NU) 

HEAD (NT) J1,J2, J3,NH, ( ( (A (I, J,K) ,1=1 ,J1) ,J=1 , J2) ,K=1,J3) 


In the forner case, the data is packed solidly into core. 
In the latter case, soae a priori knowledge of Jl, J2, and 
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J3 or their Baximum allowable values must have been 
available in order to properly dimension (A). Such a priori 
knowledge is generally contained as eleaents of (XHT) . 

Different data records may contain data of different 
types and may have differing values of J2, J3, and NH. 
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APPENDIX C 


C ARRA Y NOT ATIO N 


A shorthand notation ter referring tc arrays in the 
internal and external documentaticn of POTFAN programs has 
been developed. This notation is illustrated by the 
following examples: 


(A) This implies that A is an array. 

(A (N) ) This refers to all the words in (A) from 1 

through N. 

A (N) This refers only tc the Nth word of (A). 

(A(1,J)) This refers to all the words in the doubly 

dimensioned array A for which the first index 
varies from 1 to I and the second from 1 to 

J. 

A(X,J) This refers to the element in (A) for which 

the first index is I and the second is J. 

(A (I, J) , J=3, K) This refers to the words of (A) for which the 

first index is I and the second index varies 

from 3 to K. 

<A(I,*)) This refers to those elements of (A) for 

which the first index varies from 1 to I and 
the second index varies from 1 to some value 
which for some reason cannot be defined. 
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TABLE 3. 1.4-1 — The correspondence between IC and various 
functions defining the surface 
and default values of the functions. 


IC 

FUNCTION 

1 

XAXIS (S) 
(XAXIS (S) =S)* 

2 

YAXIS (S) 
(YAXIS=0) * 

3 

ZAXIS (S) 

• 

(ZAXIS=0) ♦ 

4 

PHI (S) 
(PHI=0) * 

5 

EX (S) 
(EX=1) ♦ 

6 

EY {S) 
(EY=0) * 

7 

EZ (S) 
(EZ=0) ♦ 

8 

YPSCAL (S) 
(YPSCAL=1) ♦ 

9 

ZPSCAL (S) 
(ZPSCAL=1)* 

10 

Not used 

IT 

V2 (V) at S = SCS(1) 
(V2=1 ) ♦ 

12 

V2 (V) at S = SCS (2) 
(V2=1) * 

• 

m 

• 

• 

• 

m 

* denotes that this is 

the default function 
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TABLE 3,1. 4-2 — Current Valid Values of lOPT and the 
Corresponding Types of Function 


lOPT 

Definition of Function V (X) 

0 

y{x) =0 

1 

Linear interpolation from the table 
(XIM(NIIi)), (YIN(NIN)) 

2 

y(x)=x 

3 

y 

4 

y=YIN ( 1) =constant 

6 

y (x) determined by controlled 
deviation interpolation method 
in subroutine CODIM. 

7 

y(x) determined by modified cubic 
spline fit in subroutine CHVFT2 

-1 

y(x) = Fourier series = YIN(1) ♦ 

YIN (2) *COS (X) ♦ YIN(3) *COS (2*x) +...+ 
XIN (2) ♦SIN (X) + XIN(3) ♦SIN (2»x) ♦ ... 

-2 

y(x) = power series in z=x-PAEAM (2) 
XIN (1) ♦XIN (2) ♦z + XIN{3) ♦z^^2 ♦ ... 
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TABLE 4, 1-1, Subprogram Call Structure 


PROGRAM 

CALLED BY 

CALLS 

ADJUST 

GEOM 

None 

AREAS 

GEOM 

None 

CFLAG 

STORGM 

None 

CHORD 1 

GEOM 

None 

CHORD2 

GEOM 

None 

CODIM 

INTRP3 

None 

CRVFT2 

INTRP3 

None 

FILEND 

STORGM 

None 

FKEAD3 

GEOM 

None 

FHOT 

GEOM 

TRANS 

BSENS1 

GEOM 

None 

FSENS2 

GEOM 

None 

GEOM 

POTGEM 

AREAS, CHORD1, CHORD2 


FREAD3, FS£NS1, FSENS2, 
NilI1, NRI2, NULLP, OBEY, 
PANL2, PBNTGM, RDFLGS, 
STSHFT, SEINEC^ SPRl, 
SPN2, SRF01, SRFIN2, 
SHFSET, STORGE, TRAN1, 
ADJUST, FRCT, 

TRAN2, XPANCP 


GRID 

PANL2 

INTRP3 

INTRP3 

GRID, SURFAS 

CODIM, CRV FT2, STRAT2 

NRI1 

GEOM 

None 

NR 12 

GEOM 

None 

NULLP 

GEOM 

None 
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CPENW 


STORGM 


None 


PANL2 

GEOM 

GRID, SURFAS 

POTGEM 

None 

GEOM, TIMEST 

PBNTGH 

GEOM 

None 

BDFLGS 

GEOM 

None 

ETSH.FT 

GEOM 

TRANS 

SETNN 

STORGM 

None 

SETNRC 

GEOM 

None 

SPN1 

GEOM 

None 

SPN2 

GEOM 

None 

SHFIN1 

GEOM 

TBAN1, TRAN2 

SRFIN2 

GEOM 

IRANI, TRAN2 

SRFSET 

GEOM 

None 

STORGM 

GEOM 

CFLAG, FILEND, 
SETNN, TIMEST, 

STRAT2 

INTRP3 

None 

SDRFAS 

PANL2 

INTRPJ, THANT, 
TRANS 

HMEST 

STORGM, POTGEM 

None 

TRANS 

RTSHFT,FROT,SURFAS 

None 

IRANI 

GEOM, SRFIN1, 
SRFIN2, SUBFAS 

None 

TRAN2 

GEOM, SRFIN1, 
SRFTN2, SURFAS 

None 

XPANCP 

GEOM 

None 


OPENW, 


TfiAN2, 
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TABLE 4.2-1 Confflon Block Usage 


COMMON 

BLOCK 

NAME 

USING 

SUBPROGRAMS 


C0M1 

GEOM, 

GRID, PANL2 


COM2 

GEOM, 

PANL2 


COM3 

6E0H, 

PRNTGM, STORGM 


COM4 

GEOM, 

RTSHFT 


CONST 

POTGEM 

FILEND 

RDFLGS 

SURFAS 

, GEOM, CFLAG, CHORD1, CHORi)2, CHVFT2, 
, FREAD3, GRID, INTRP3, CPENW, PRNTGtt, 
, SPNl, SPN2, SRFIN1, SRF1N2, STORGM, 

, TRAN1, TRAN2, XPANCP, ADJUST, FHOT 

SREDAT 

GEOH, 

SRFIN1, SRFIN2, 

SRFSET, SURFAS 
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TABLB 4.3-1 Logical Units Used by PCIGEM 


FORTRAN 

VARIABLE 

LOGICAL UNIT DESCRIPTION 

DEFAULT 

VALUE 

NTCP 

Output device for 
conversational prompts 
and error messages 

b 

NTP 

Output device for 
normal printout 

6 

NTR 

Card input device 

5 

NTS 

Device for storing geometry 
file 

1 
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TABLE 7.7-1 Values of the Mean Camber line 
of the NACA 5-Digit 230 Airfoil 


X 

Z/c 

-.0200 

-.0067692 

0.0000 

0.0000000 

.0200 

.0054767 

.0400 

.0097886 

.0600 

.0130632 

.0800 

.0154283 

. 1000 

.0170115 

. 1200 

.0179405 

.1400 

.0183864 

.1500 

.0183864 

.1600 

,0183463 

.1800 

.0189785 

.2000 

.0176670 

.2025 

.0176119 

1.0000 

0.0000000 
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FIGURE 3.1-1. Data Needed to Define an. Arbitrary Surface. 
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FIGURE 3.1.6-i. Method of Calculating Unit Normals 
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0 = CORNER POINT 
X = CONTROL POINT 

1 I I 
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FIGURE 3. 2. 2-2. Continued. 
Figures-8 
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FIGURE 


.2,2-2. Concluded. 
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V, K2 


SUV) 


SU(V) 


FIGURE 3. 



.3-1. Definition of Various Quantities Required in the Calculation 
of (S,V) Values in the Interior of a Segment. 
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XGPS 


XGPVU(Kl) 



XGPSU(K2) 


FIGURE 3. 2. 3-2. Determination of FS(K1,K2) and FV(K1,K2) from the Non- 
dimensional Description of the Grid Line and Segment Boundary 
Intersections. 
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FIGURE 3. 2.3-3. Example of Nonuniform Network of Corner Point 
(Connected by Line Segments) and Boundary 
Condition Point (shown as dots) Grid Lines. 
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FIGURE 4.1-1. POTGEM Program Subroutine Structure. 
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DSFLAG 


VORTEX MODEL 


DSFLAG 


VORTEX MODEL 
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FIGURE 6.1-1. Vortex Models Available in VVIM (Cont'd). 
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DSFLAG 


VORTEX MODEL 


DSFLAG 


VORTEX MODEL 



FIGURE 6.1-1. Vortex Models Available in VVIM (Cont'd). 
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DSFLAG 


VORTEX MODEL 


DSFLAG 


VORTEX MODEL 



FIGURE 6.1-1. Vortex Models Available in VVIM (Gont'd) 
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Y, S, N1 

■3° 



X,V, N2 (1.7071,1.) 

ASPECT RATIO = 2sj2 
TAPER ratio = 1/3 

FIGURE 7.1-1, Thin, Symmetrical, Swept, Flat Wing (Warren 12 Planform) 
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1, T 

2, TITlf 

3, tfst case 1 ; Warren u winc 

fl, WlNp 

5, ^>OATA CROOT«1.0S07,rTlP«,5536,L*MLe«55,5fl3,Pa»l 

6, nsEGMENTS 

7, tOaTA NflPS»lo,N8PV«a SEND 

VL8C 

R, sOaTA send 

lo; Slac 

It, tOaTA I0PT»? stNO 

!2, PaN| 

13, tOaTa RSlsT SEND 

l«, osel 

15. I .1 

16. 1 .1 

1 7’, 3G 

18.0 

IR, finish 

20, tOATA LOGfl2)*t,T,lNTfl0)«l,ELT(5)»1.0,0,9,0,0 

21, store 

22, ..OaTA lOal SEND 

23, PWIMT 

2«, -tOATA PR1 nT*1s*T sEnD 

25, STOP 


FIGURE 7.1-2. Input for POTGEM Test Case 1. 


tO send 


send 
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POTFan tEUMfeTRY PRoGRAw. VERSTON 1 3 

ssvssBaoaassMsaaaaaaaasscasPawsPaaas 

TIME m 08/0R/76 02«23JS« 

ENTER BATCH 
♦TITIF 

test case 1 . WaRREN 12 wING 


+ WlNr, 

♦DSECMENTs 
+ Vl RC 

♦ SLBC 

♦ PANl 

♦dsfl 

♦FINISH 

♦store 

FiLg i.cM-PNC/Lie* has Been opeNeo fdr writing on unit \ 

CREATION TIME * 08/0R/76 02 j 2<II0S 

creation uf Geometry fiuf 


title « TE^^T case I . ,,aRREN 12 ^ING 

CLOG) « fffftffffffttffftftf 

(INT) » oil Bio ^ 0 0 0 

CELT) s .70715000 r.OOOOOOO '70715000 

iSgOfelUSI-.O? 0 . o', 

♦print 

PRINTOUT UF Geometry ftIe data 


TITlE * TEST CAgE 1 . WARREN 12 WINr 

creation time a OB/OR/76 02l20io5 
(IFORM) a noillllll 


FIGURE 7.1-3. Output for POTGEM Test Case 1 


1 0 
}' 0000000 

0 . 


1,0000000 
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no) 

■ 

1 







(lOG) 

■ E 

F F F T F r F 

F F F T T F F 

F T F T F 




w 

(INT) 

s 

Oil S 

lo 0 0 

0 0 

0 

loo 

0 

CfLT) 

• 

.7071300 

1. 0000000 

7071500 

r.OOOOOOO 

r oooooo 

8| 



1 3'. 298 1080 

'0013298 

0. 


o'. 

0. 










_ K 

PiNEl 

C0«N-f 

H points 







I 

J 

XfT,J) 

Y (I * J 1 

7rT»J) 


Sfl.Jl 

V(T,J) 

»E P> 

1 

1 

;o882937 

•.0000000 



o’. 

.0862938 


? 

1 

M972289 

'1000000 

o' 


lOOOOOO 

.1972289 

W ?ri 

1 

1 

'3281801 

,2000000 

o' 


’.2000000 

.3201891 

^ Ft 

a 


'0500993 

3000000 

0' 


' 3000000 

.9590993 


5 

1 

'5900305 

'aoooooo 

o' 


.9000000 

.5900395 


6 

1 

'7209897 

'soooooo 

o; 


’.5000000 

.7209897 


7 

1 

'8519009 

,8000000 



>000000 

'.0519099 


B 

1 

'9820001 

,7000000 



7000000 

.9828901 


q 

1 

l' 1137753 

,8000000 

o' 


>000000 

1.U37753 


10 


l'200710« 

,9000000 

o' 


>000000 

1.2997105 


11 

1 

l' 5758057 

i'ooooooo 

o; 


1.0000000 

r.3758957 


1 

? 

'3310680 

..0000000 

o; 


o'. 

'.3319668 


? 

2 

'0007280 

' loooooo 

0 . 


'.lOOOOOO 

.9997289 


3 

2 

'5579801 

,2000000 

o; 


'.2000000 

.5579891 


a 

2 

'8712018 

,3000000 

o'. 


>000000 

.8712910 


5> 

2 

^7800995 

,9000000 

0 ' 


>000000 

. 7899995 


6 

2 

'8977572 

.5000000 

o; 


’.5000000 

.8977572 


7 

2 

I'oiiOioq 

'8000000 

o'. 


.8000000 

1.0110199 


8 

2 

!' 1202728 

'7000000 

o'. 


'. 7000000 

r. 1292728 


q 

2 

l'2375303 

,8000000 

o’. 


'.8000000 

1 .2375303 


10 

2 

l'3507880 

.9000000 

o’. 


>000000 

1 .3507880 


li 

2 

r 0800057 

r.OOOOOOO 

o; 


r.OOOOOOO 

1.9890957 


1 

5 

'5988037 

..0000000 

o'. 


0. 

.5988938 


2 

3 

'8922239 

' lOOOOOO 

o'. 


'. lOOOOOO 

'.8922239 


3 

3 

'7878001 

'2000000 

0 ' 


'.2000000 

,7878091 


4i 

3 

'8833803 

,3000000 



'. 3000000 

.8833893 


5 

5 

'9789805 

.9000000 

o'. 


.9000000 

,970<>895 


6 

3 

r 0705007 

'5000000 

o; 


.5000000 

1 .0795997 


7 

5 

t'l70l209 

,8000000 



>000000 

1.1701299 


8 

5 

1'2857o51 

,7000000 

K 


' 7000000 

1.2857051 


q 

3 

1 ' 3812853 

,8000000 

0 . 


>000000 

1.3812853 


10 

3 

i'058885O 

.9000000 

o’. 


'.9000000 

1.9588855 




FIGURE 7. 

1-3. Output 

for POTGEM 

Test 

Case 1 (Cent 

d). 
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11 

3 

t'5524457 

l' 0000000 

O’, 

I’.OOOOOOO 

1.5524457 

1 

a 

'8S18J87 

.,0000000 

o: 

0. 

.8618188 

2 

a 

'<J5R7214 

MOOOOOO 

O'. 

1000000 

.9597214 

3 

4 

r 0178241 

^2000000 

O’. 

’.2000000 

1.0176241 

a 

tt 

1^09S52t»8 

3000000 


; 3000000 

1 .0955268 

5 

4 

1^1754295 

^4000000 

o: 

’.4000000 

1.1734295 

6 

u 

1'2S13322 

,5000000 

o'. 

'5000000 

1.2513322 

7 

4 

1^3292349 

6000000 

o' 

.6000000 

1.3292349 

B 

4 

1'*4071376 

'7000000 

O’. 

'7000000 

1 .4071376 

<J 

4 

1'4850401 

,8000000 

0. 

.8000000 

1 .4850403 

10 

4 

r'Sfc29430 

,9000000 

O’. 

'9000000 

1 .5629450 

11 

4 

t'6408457 

1. 0000000 

o; 

1.0000000 

t .6408457 

1 

5 

1^0607000 

..0000000 

o: 

0. 

1.0607000 

2 

5 

r 1253445 

'1000000 

0. 

'.1000000 

1.1253446 

3 

5 

l' 1899891 

,2000000 

o; 

'.2000000 

1.1899891 


S 

r 2546337 

,3000000 

0 

/5000000 

1,2546337 

5 

5 

i'3i92783 

,4000000 

0 

.4000000 

1.3192783 

6 

5 

1.5839228 

,5000000 

o'. 

’.5000000 

1.3639228 

7 

5 

r 4485674 

,6000000 

o'. 

.6000000 

1,4485674 

a 

5 

l'5l32l20 

,7000000 

O'. 

'7000000 

1.5132120 

9 

5 

i'5778565 

,8000000 

o; 

.8000000 

1.5778565 

10 

5 

r 6425011 

9000000 

o'. 

'9000000 

1.6425011 

It 

5 

1^7071457 

1.0000000 

o’. 

1.0000000 

1.7071457 

Unit 

veltors along Wake 

ELEMENTS 




I 

J 

tiVwX(T, J) 

UVWYfl.J) 

Ijvw7nf J) 




unavailable 


boundary condition el AGS 
I J BCFlAKfl.J) 

unavauabie 


doublet singularity ELaGS 
I J DSFLAG(I#J) 

1 1 5o 

2 I 5o 


FIGURE 7.1-3. Output for POTGEM Test Case 1 (Cont'd). 
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FIGURE 7.1-3. 


1 1 
« 1 

5 I 
*> \ 
1 1 

8 I 

9 I 
10 1 

t 2 

2 2 

3 2 

4 2 
« 2 

6 2 

7 2 

8 2 

9 2 

10 2 

1 3 

2 3 

3 3 

0 3 

5 3 
t> 3 

7 3 

8 3 

9 3 
io 3 

1 4 

2 « 

3 « 

4 4 

5 4 

6 4 

7 4 

8 4 
4 

10 4 


3o 

5o 

30 

3o 

3n 

3n 

3p 

3o 

3n 

3n 

3o 

30 

3o 

3o 

3o 

3fi 

5o 

3n 

3n 

3ft 

3ft 

3ft 

3ft 

3ft 

3ft 

5fl 

3ft 

5ft 

3ft 

5ft 

3ft 

3ft 

3n 

3ft 

3ft 

3ft 

3ft 

3ft 


SOURCE SINGULARITY FLaGS 

Output for POTGEM Test Case 1 (Cont'd). 
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I J SSFtAGfl.J) 


UNAVAIUABIE 

boundary condition poinis 

I J .Hr i-i. J) 

VrC(I»J) 

ZBCf I. J) 

0 ; 

SbC{I,J) 

VRCf I,J1 

1 

2 


'259929S 

' 0500000 

’osooooo 

.25RR25S 

t 

' 3820259 

1500000 

0 ' 

,1500000 

.3820259 


1 

^5091223 

'2500000 

0 : 
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FIGURE 7.1-3. Output for POTGEM Test Case 1 (Cont'd). 
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FIGURE 7.1-3. Output for POTGEM Test Case 1 (Cont’d). 
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FIGURE 7.1-3. Output for POTGEM Test Case 1 (Concluded) 
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FIGURE 7.1-4. Planview of POTGEM Test Case 1. 
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FIGURE 7.1-5. PLOTGM Input That Generated Figure 7.1-4. 


Figures-33 



T 

T!T|f 

TfST CftSE ? ; WaR«FN 12 wimr 

C«»r 

Swy 1 

•-ImCRVI IC = 1,COPT = 0 $En» 

St' II 

UncRVI IC*?.rOPTs? iEND 

Si-'T! 

lCS«,COPTsi,NTAfl*I, VARgSK'-ifK, $ENd 

st^n 

ilNCRVl IC=S,COPTcO $EnO 
SRI! 

tINCRVJ ICs7,C0PIsl , VA«2*t. $FNn 
SflT j 

tI^,CRVl ICsil.COPIsO ttNf) 
nsEGMfNT.S 

-GaTA NBPfj*to,NBPV*« ?ENO 
VLBC 

^OaTA *FND 
SLPC 

iDATA lOPT*? «ENd 

sl 

-Data inPTSv*n hEnD 

sU 

s-DaTA IOPTSV*j ,NTABSV»i ,VAR2SV«1 _0 *ENd 
VL 

•:DaTA NtABSV»2,VAr1sV*0.,1,,VAr 2SV*0,, 1.3535*157 jENn 

vu. 

iDaTA VAR?SV*1 ,06 oT, 1 ,7o7ia57 $E;ND 
PA.Nl 

,DaTa RS1*T $EN() 

D'^Pl 

i .1 
1 .1 
so 


FIGURE 7.2-1. Input for POTGEM Test Case 2. 
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FIGURE 7.2-1. Input for POTGEM Test Case 2 (Concluded). 
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FIGURE 7.2-2. Output for POTGEM Test Case 2 (Cont’d). 
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FIGURE 7.2-2. Output for POTGEM Test Case 2 (Cont’d). 
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FIGURE 7.2-2. Output for POTGEM Test Case 2 (Concluded). 
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FIGURE 7.3-2. Input for POTGEM Test Case 3 (Concluded). 
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FIGURE 7.3-3. Output for POTGEM Test Case 3. 
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FIGURE 7.3-3. Output for POTGEM Test Case 3 (Cont’d). 
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Output for POTGEM Test Case 3 (Cont'd). 
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FIGURE 7.3-3. Output for POTGEM Test Case 3 (Cont'd). 
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FIGURE 7.3-3. Output for POTGEM Test Case 3 (Cont'd). 
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20,5200000 

3 

16 


o'. 

.1,0000000 

20.3200010 

nTOP VH Tr.Rh 




1 

J 

NtOPXfl, J) 

NTOPYf l»>n 

NT0P7{T, J) 


unavailable 




NBOt vFt-TnRjj 




I 

J 

N, OTXf l.J) 

NBOTVfl.J) 

NROTZf I,J) 



unavailable 


VELOCITV *LONg NTOP VEcTORS 

FIGURE 7.3-3. Output for POTGEM Test Case 3 (Cont’d). 
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! J 


«T0P{!» J) 


UNAVAlLABir 

VELOCITY »tONr. NBOT VECTORS 
I J OHOT(l,J) 

unavailable 




CORNER ROINTS ALONG VL ANU VU EDGES 


I 

XvLrf 1) 

TVLCM) 

?VLcf i> 

SVLCTIJ 

VVLcd) 

1 

.10.00000 

0. 

.6,09600 

;iO, OOOOO 

.6,09600 

2 

•,00000 

0. 

•6.09600 

•,00000 

•6,09600 

S 

10' 00000 

0, 

.6,09600 

10,00000 

.6,09600 

a 

20' OOOOO 

0. 

.6.09600 

20,00000 

.6,09600 

bounoary Points 

ALnNG VL and 

VU EDGES 



I 

XVLm(I) 

VVLBfn 7VLB(n SVLBCD VVLbCU 

xvUBf n 

1 

•5,000 

0 . -0 

,096 .S, 

000 .6,096 

.5,000 

2 

5' 000 

0. .6 

.096 5, 

000 •6,096 

5,000 

3 

15'000 

0. .8 

,096 15, 

000 .6,096 

15,000 

roRNtfi 

t points along SL *nD su 

EDGES 




XvSLf (I) 

YSLCfl) 

ZSLCfl) 

SSLcfl) 

VSLc(l) 

i 

.10. OOOOO 

0. 

•6,09600 

;iO, OOOOO 

.6,09600 

?. 

•10.00000 

2,03200 

.6,09600 

.10,00000 

.0,06000 

3 

•10.00000 

0,00400 

.6,09600 

^10,00000 

.2,03200 

u 

•10,00000 

O.O’OOO 

.6,09600 

.10,00000 

.90,00000 

s 

.10.00000 

7,07977 

.5,79760 

^10,00000 

^72,00000 

6 

•10,00000 

9,67910 

.0,95177 

•10,00000 

;S0, OOOOO 

7 

•iOtOoOOO 

11.02777 

-5,583l0 

•10,00000 

.36,00000 

8 

•lo.ooooo 

11.89560 

-1,88577 

.10,00000 

.18,00000 

9 

-10.00000 

12,19200 

^,00000 

•10,00000 

.,00000 

10 

-10,00000 

11 89360 

1,88377 

.10,00000 

18' OOOOO 

11 

.10,00000 

11^02777 

3.58310 

.10,00000 

36.00000 

12 

•10,00000 

9,67910 

0,95177 

.10,00000 

50.00000 

13 

•10,00000 

7 97977 

5,79760 

.10,00000 

72.00000 


xvuc'rll 

wucd) 

zvucflJ 

svucvl) 

vvurri) 

.10,00000 

o: 

6,09600 

.10,00000 

.6,09600 

.,00000 

O’, 

6,09600 

.,00000 

.6,09600 

10’ 00000 

0. 

6,09600 

10,00000 

,6,09600 

20.00000 

K 

6,09600 

20,00000 

•6,09600 


YVtiBf n 

zvuBf n 

SVUBin 

VVURf I) 

C 0 «D 2 f n 

SPaNIcH 

o; 

6,096 

.5,000 

-6,096 

0. 

0. 

0 

6,096 

5,000 

.6,096 

0. 

0. 

o: 

6,096 

15,000 

.6,096 

0. 

0. 


xsucf n 

YSUcfl) 

ZSUcfl) 

ssucf n 

vsuc(i) 

ao'.ooooo 

0. 

.6,09600 

20,00000 

.6,09600 

20’ OOOOO 

2,05200 

•6,09600 

20,00000 

•4,06000 

20,00000 

0,06000 

•6,09600 

20,00000 

.2,03200 

20,00000 

6,09600 

•6,09600 

20.00000 

.90,00000 

20^00000 

7,97977 

-5,79760 

20,00000 

.72,00000 

20,00000 

9,67910 

-0,93177 

20,00000 

.50.00000 

20,00000 

11,02777 

-3,58310 

20,00000 

.36,00000 

20,00000 

11,89560 

-1,80377 

20,00000 

.18,00000 

20,00000 

12,19200 

.,00000 

20,00000 

•,00000 

20,00000 

ir. 89560 

1.80377 

20,00000 

18,00000 

20.00000 

11.02777 

3,50310 

20,00000 

36,00000 

20’ OOOOO 

9.67910 

0.93177 

20.00000 

50,00000 

20.00000 

7,97977 

5,79760 

20,00000 

72,00000 


FIGURE 7.3-3. Output for POTGEM Test Case 3 (Cont'd) 
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« 1 0 • 00000 

6.00600 

6,09600 

.10,00000 

o ’. 

20' 00000 

6.09600 

6.09600 

20.00000 

0. 


.10.00000 

4.06400 

6,09600 

.10,00000 

.2,03200 

20.00000 

4.06400 

6.09600 

20.00000 

.2.03200 

16 

.10,00000 

2.03200 

6,09600 

.10.00000 

.4,06400 

20.00000 

2.03200 

6.09600 

20^00000 

.4,06400 

17 

.10,00000 

0. 

6,0o600 

.10,00000 

• 6,09600 

20.00000 

0. 

6,09600 

20.00000 

.6,09600 


BoUnOaR'^ ROInTS aLonG 8L AfjD SU EDGES 


I 


VSLBfll 

? 8 LBf n 

SSLBfil 

VSLB(l) 

XSUb(I) 

T8Ub(1> 

7 SUB ( 1 ) 

SSUb(U 

VSURfl ) 

CORDl f I ) 

SPANEfU 

1 

. 10.000 

1.016 

. 6.096 

. 10.000 

. 5,080 

20.000 

1’ 016 

— 6 , 096 

20.000 

. 5,080 

0. 

0. 

2 

. 10,000 

3,048 

.6 , 096 

. 10.000 

. 5,048 

20.000 

3.048 

. 6,096 

20,000 

. 3,048 

0. 

0. 

3 

. 10.000 

® i , 0«0 

. 6,096 

- 10,000 

. 1,016 

20.000 

5 ‘ 080 

• 6,096 

20,000 

-1,016 

0. 

0. 

4 

. 10' 000 

7.050 

. 6.021 

-lOlOOO 

. 81 '.000 

20.000 

7,050 

— 6 , 02 1 

20,000 

. 81,000 

0. 

0 . 

*) 

. 10.000 

0,864 

. 5,432 

. 10,000 

. 63,000 

20,000 

8.864 

- 5,432 

20,000 

• 63,000 

0. 

0. 

6 

. 10.000 

10.407 

. 4,311 

. 10,000 

- 45,000 

20.000 

10 ’! 407 

- 4,311 

20,000 

- 45,000 

0, 

«. 

7 

. 10.000 

11,528 

. 2,768 

. 10.000 

- 27,000 

20.000 

11,528 

. 2,768 

20,000 

- 27,000 

0. 

0 . 

8 

. 10.000 

12.117 

.,954 

. 10,000 

- 9,000 

20.000 

12 ; il 7 

;,954 

20,000 

- 9,000 

0. 

n. 

9 

. 10.000 

12.117 

,954 

. 10,000 

9.000 

20.000 

121117 

,954 

20.000 

9.000 

0. 

0 . 

10 

. 10,000 

11 '.528 

2,768 

. 10,000 

27 000 

20.000 

11,528 

2,768 

20,000 

27.000 

0. 

0 . 

11 

. 10,000 

10.407 

4,311 

• 10,000 

45.000 

2o;ooo 

10.407 

4,511 

20,000 

45.000 

0. 

0 . 

12 

. 10,000 

8.864 

5,432 

. 10,000 

65.000 

20.000 

8 ' 864 

5.432 

20,000 

63,000 

0. 

0 . 

13 

. 10,000 

7 050 

6,021 

. 10,000 

81.000 

20.000 

7'050 

6,021 

20.000 

81.000 

0. 

0 . 

14 

. 10.000 

5,080 

6,096 

. 10,000 

. 1,016 

20.000 

51080 

6.096 

20,000 

- 1 ,016 

0. 

0 . 

IS 

. 10.000 

3,046 

6,096 

. 10,000 

. 3,048 

20.000 

31048 

6,096 

20,000 

. 3,048 

0. 

0 . 

16 

. 10.000 

!>16 

6,096 

. 10,000 

. 5,080 

20.000 

r . oi 6 

6,096 

20,000 

. 5,080 

0. 

0 . 


Force StRSlNG tOCATIONS in nJ»01RECTIon 
I J *SlfT»J) VSJfTfJ) 

UNAVAIIArU 


FORCE sensing locations IN N2-DIREcTlON 

T J *vS?fTtJ> ^SH(TfJ) 2S2fl»J) 

UNAVAILABLE 
♦ STOP 
STqP 777 

FIGURE 7.3-3. Output for POTEGM Test Case 3 (Concluded). 
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NRBR RME5 fOXSO TUNNEL 

FIGURE 7. 3-4 (a). Top View of POTGEM Test Case 3. 

Figures-59 



Figures- 60 


FIGURE 


7.3-4(c). Oblique View of POTGEM Test Ease 3, 
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TRUE 

THIS IS A PLOTGM RUN TO PLOT AN X-Y VIEW/ Y- Z VIEW/ AND 
OBLIOUE VIEW OF A SAMPLE PANELLING ARRANGEMENT ON THE 
NASA AMES 12.2 M. BY 24.2 M. (40X30) WIND TUNNEL. 

THE PANEL ARRANGEMENT TO BE PLOTTED IS THAT PPODUCEP BY 
POTGEM TEST CASE 3 AS DESCRIBED IN NASA TM-X 73/ 127/ 

"NASA AMES POTENTIAL FLOW ANALYSIS (POTGEM) GEOMETRY 
PROGRAM (POTGEM). VERSION 1." 

READ IN THE GEOMETRY FILE 

READ 

3 

INITIALIZE THE PLOTTER FOP. AN OFFLINE PLOT AT 15 
CHARACTERS PER SECOND. 

THE PLOTTER DIRECTIVES WILL BE STORED ON FILE 3 . PT FOP LATEP 
PLOTTING. THE FILE WILL BE CREATED BY WRITES TO UNIT 3. 

IPLOT 

15 3 8 

PLOT THE PLAN VIEW (X-Y VIEW) 10 CM. TO THE 
RIGHT AND 10 CM. ABOVE THE CURRENT PEN LOCATION. 

PLOT 

$DATA YOFF=10/XOFF=10/XSCALE=1 .0/YSCALE=l .0 $END 
NOW PLOT THE CONTROL POINTS. THESE WILL 
SHOW UP AS DOTS. 

PBCP 


FIGURE 7.3-5. PLOTGM Input That Generated Figure 7.3-4. 

Figures-62 


WA KE S 

$DATA WAKECf^=5 .0 ,A RRWC M= . 75 $END 
ENTER T ITLE AND PLOT IT. 

RT ITLE 

NASA AMES 40X80 TUNNEL 
PTITLE 

$PATA XT ITLE=10,YT ITLE=5,HE IGHT=. 762 SEND 
NOW PLOT THE FRONT VIEW (Y-Z VIEW). 

PAGE 

rzviEw 

PLOT 

$DATA $END 
PBCP 

FINALLY, PLOT AN OBLIQUE VIEW OF THE TUNNEL. 

PAGE 

OBVIEW 

PLOT 

$DATA $END 
PBCP 
WAKES 
$DATA $ENO 
STOP 


FIGURE 7.3-5. PLOTGM Input That Generated Figure 7.3-4 (Concluded). 

Figures- 6 3 











2. 

«. 

5 . 

6 , 

7. 

8. 

10 . 

n. 

12 . 

13 . 

15. 

16 . 

17 . 

18. 


rv, 

21 . 

22 . 

2i. 

2«. 

25 . 


26 . 

27 , 

28, 

30 . 

31 . 

32 . 

”, 

3«. 

35 . 

36. 


ItTl? 

TfST CASE a i,- IHIN, SWEPT, UNCaHBERED, UNTWISTED WING WITH DlHEORiL 

C»RY 

SPTJ 

»TNCRVl 1C»1,C0PT*3 $END 
Su- T 1 

tlNCRVl 1C»1 I ,COPT»0,SCSsO. $ENO 
SPTi 

-iINCRVI iCB!2,NTAfiBj,VAR2f 1 )«0,5,Cf)PT*l,SCS«2 .iEnD 
nsEGMENTS 

-.0*T* NrP 3(. )s2.NRPV(, )»2 send 
VI HC 

«oata send 

S^RC - 

,OaTA IOPTs? sEHO 
®L 

•«OaTA IOPTSVbo send 
RU. 

<i.OaTA IOPTSVbI ,NTABSV*1,V*R2SV( 1 )«2 $EnO 
VL 

tDaTa NTAHSV»2,VARlSVf i)»o,2,VAR2SV(i)*o,o,5 SEN(> 

VIA, 

fOAT* VAR2SVn)*2,1.5 sEwD 
GRlp 
PaN| 

+ OaTa RSl*T,RS2*T sENf) 

RaSS 

*DaTA EAXiSn 1*0,0 , -I.PHIsRO TEND 
RrSS 

vDaTA EAXIs( lleO, 1 ,0,PHl«a5 *EmD 
DsFl 

1 2 
1 

2 

I 


FIGURE 7.4“4. Input for POTGEM Test Case 4. 
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37. 

? 


3B. 

h 


39. 

?. 


«o. 

? 



12 


«2. 

0 


«3. 

rlMlSH 


fl«. 

»data 

tOC(lS)«T, 1NT(10)*1,FLT(5)»1»0,0,FLT(I )»5 sEMP 

«5. 

STORF 


<16. 

rDATA 

10»a 5tND 

«7. 

print 


«e. 

iOaTA 

PRINT{l)tlS*T SEND 

«9. 

STOP 


FIGURE 

7.4-4. 

Input for POTGEM Test Case 4 (Concluded) . 


Flgures-68 




POTFAN UfeiiMMPV PROGRAM, VtRSIUN l'.? 

«ssi:aas»sssz«cS3«s«ss3sa»s«saK8S3s£ss 

T1*^L s 0H/ofV7<’ 

eNTH- RATt H 
M I T! F 

TtST CA>,F 'J -- T^'IN, ijWF.PT, UNCftMHERFO, UfJTWlsTtP “UNG W!t^ OIHEORaL 


»CAKV 

♦ SR 1 t 

♦ 5H 1 1 

♦ SRM 

+ OvStr,HFNT* 

♦ Vl.hf 

♦ SlPf 

♦ 5L 

♦ S'l 

♦ Sfl 

♦ Vll 

♦ GRIP 

NStGS s < NSFT.V s I 


XGPVt (. * 

-1,0000 f>. r.oooo 


XGPVI'C. t 

"1.0000 r.oooo 

XGPvU^ -*• 

".bOOO ,'-000 


XC,Pvllfi r 

-.boon ,M)oo 


MbPSfNSfcGSi s P NRPVfNSEGV) r 


? 


FIGURE 7.4-5. Output for POTGEM Test Case 4. 
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XGPSI r * 


7S00 

, <•* S A 0 

l'. 0000 

xr.pstic f 



-,7'^00 

■.-’500 

r.oooo 

XGPhI R = 



-,?S00 

. /500 


xGPSHR * 
-.<JS00 

fRArji 

tWOSS 

+ ROSS 
♦ pSM 
tPlNlSM 
+ STi;S{. 

. ?5oO 



FiLf r,/l. IH^ HAS been OPeNKpi FUR WHITING ON UNIT i 

CHFaTIos. »i-€ t 07»5a»^»4 

UN11 5 (•-’•n'lLED AND R^lEAStn 

CRFaIION i>F nfPHFTRV MLT 


TlItF; « •F^T f*Sf <i ThIN, SWEPI, UNCamHEREO, IJnTwISTEd wIng wITh nlMtOWAl 


(LOG) s 
( I H n = 

t f- F F T F F F F F F r T F F F 
0 r \ i 2 0 

T F 1 1 
0 0 

n 1 0 

no 




fFl T 1 = 

♦ PR 1 *»T 

T.OQOOOOO ^.OOOOOOO 1 

.'*6^025a0F-0? 0, 

1 '.0R06601 
O', 

i' 0000000 

0. 

1.0000000 

0. 

o'. 

Sfe, 6025^0 


PHTnIUUI it G^OMETpy FILE LUTa 


TITiF TEST CASE '1 — ThIm, SWEPT, UNC a-^HEHED, UnTRISTEO wUjG WITH PIHEDHaI 

CrEMIUN tl*-F S OB/09/76 ottiatPa 
f irCiRH) = 11011111111 
flP) = a 

EIGURE 7,4-5. Output for POTGEM Test Case 4 (Cont'd). 
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ftOfr) 

s F f- 

F F T F 2 F 

F F F F T f F F 

T F . t T 



UNT) 

z 

0 5 5 

2 2 0 

0 0 0 

1 0 0 

ft 

(FLt ) 

» 

3.0000000 

a'.ooooooo 

! .0606601 

2,0000000 

1 ,0000000 



«6.6021>000 

;oo«*&05 

0. 

0. 

0. 

PAWfi 

CUfiNf-t- 

fo|nT.s 





I 

J 

X f I * J 1 

y f I , J) 

?cT»vn 

sn # Ji 

V n , J ) 

1 

1 

; ! 7b7767 

..0000000 

«i, l 7677^7 

o’. 

,2500000 

2 

1 

' ISO 

I'.OOOOOOO 

-, 15?58?5 

1 .0000000 

,0375000 

3 

1 

' 7'?S«'»5l 

2.0000000 

,,0883863 

2.0000000 

.62S0000 

1 

? 

'8238855 

.,0000000 

,.8838835 

0. 

1.2500000 

2 

2 

|'^0»««<»F50 

r.OOOOOOO 

-, 6620126 

I’.OOOOOOO 

1 .1875000 

3 

2 

1 ' 1«‘»0085 

2'0000000 

-.aai'juiF 

2.0000000 

1.1250000 

1 

3 

r‘«l02l5S 

.,0000000 

,1 ,aia2i 55 

o’. 

2,0000000 

2 

3 

roia2t3S 

1’ 0000000 

■ 1 , 06Q6602 

r.OOOOOOO 

1.7500000 

i 

3 

)'0V02i 35 

2.0000000 

., 7071068 

2.0000000 

1.5000000 

UN! T 

VECTiiRri 

ALO^'I. KakF 

EtFMFMTS 




I 

J 

Vvtxn, j> 

hVWYM , J1 

iiVF2( Tf J5 




UNAVAH Abl fc 


B 0 UNt 5 A»* i.fiWOnrnN KARS 
I J HC FI ARC l»s>) 

unavaufhi I 


OOUhUT SJNC.UI sBITv fJ aOS 

1 J PSFLAR(I,J 1 

11 ^ 

2 t ? 

12 ^ 

2 2 12 


O C> 

is 

II 

pi ^ 



FIGURE 7.4-5. Output for POTGEM Test Case 4 (Cont'd). 
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SOURCfc i»l TV FLAOti 

1 S F (. A fi f I » J ) 

UNAVAUFBIF 


BOwNOaRV 
I J 

1 0‘'i>TT!ON BQlwtS 

VhC(TfJ) 

/8Cf I,J) 

ShC ( T, J) 

VBCfl.J) 

1 

1 

'hryoBiBS 

Isoooooo 


’.soooooo 

,781£500 

? 

1 

''H6J 786U 

r.5000000 

« , 33 ! <*5t> 3 

J',5000000 

.H4J7500 

1 

? 

r‘«?S9S$39 

' soooooo 

,1 ,0»27ST2 

.5000000 

1 ,6562500 


e 

i'i0i7?81 

1.5000000 

7753980 

1 .5000000 

1 . 4687500 

liNH 

I 

normal b and areas 
J '<NxfT«JT 

HNV f T ,.n 

tlN7 r I . J> 

Dat I.JI 


1 

1 


Ua^S■^Sb 

■'» , 6BS99o 

.9019^94 


? 

1 

-,E>RS99« 5 

20?S3Sb 

-■ ..b859<»<15 

.9U4?35i 


1 

e 


.2«2'jS56 

,,b»599il5 

, 9764470 


? 

? 



.,08^9941 i 

. 4851 764 


MTOP 

T 

VEf T 
J 

ik.'s 

NtnpXf 1, J) 

NTOPVf t. J) 

NtOPV( T, J) 




U(^AVA iLAm t 


N«0! VFCT'R.- 

I J N. nixfi,j) Nuorvfi,j> v.,0T/fi»J) 

liNavailabie 


VELOrn^ f:LONo NTOP VfccTORS 
I J f'TriPfi,.!) 

URAVA U ABt.E 

FIGURE 7.4-5. Output for POTGEM Test Case 4 (Cont’d). 
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VFLOrtT^ NOOT VfcCinRS 

I j unnTf J f ,J ) 

UNAVAILAHtf 


rO«NfK >»oiNTS &t oNr, vi ano vu fdrks 


I 

XvLrfTl 

vvLr f I) 

7VLrf n 

SVLc f n 

vvi r f t ) 

1 

0' 

0. 

o' 

0. 

o'. 

2 

’ 5*'5S'4 

1,00000 

.00000 

1 ,00000 

.25000 

5 

' 707 t 1 

2.OO000 

'.00000 

2'.00000 

]soooo 

HOilNtH«Y t pints 

along V(_ a no 

vu EOGfcS 



I 

KVLr:( 1 ) 

¥vi Mf n £vi nf n svi-pf 1 1 

VVLHf I ) 

xviirifn 

1 

' r/7 

.500 

.000 .-’-500 

.12S 

1 414 

2 

’ 5 <0 

l.SOO 0 


’. 5/5 

1 |4iu 

COHnf H 

points alPXic: SL ANf) SU 

FDGPS 



T 

X j f J ) 

¥S| cf n 

7S| rrl> 

5S| rf n 

vsi r f n 

1 

1 ?f> 7 « 

-.00000 

17P78 

0. 

'.25000 

2 

' 8}*5«8 

-.00000 

-.86588 

0. 

1 *25000 

4 

1 ' a 1 4 ? 1 

..ooooo 

-1.4142! 

0. 

2.00000 

«nuNtuwi^ Points 

*lqng si and 

SD EOGtS 



T 

XSl. ( I ) 

YSlUfM /SlMfl) SSLP.fl) 

VSLRf 1 1 

xsunf n 

1 

.S‘-0 

-.000 

,550 0 

.7S0 

9 72 

2 

r 2 47 

-.000 -1 

. 25 ? o' 

l'.7.;rt 

1 . 52 h 


xvlirrn vv.irrn sviicfl) wncfl) 

r.£(l«?l -,00000 0, 00000 

I'uju^l l-jOOOOO wt.OeOtoo 1,00000 1.7S0O0 

r’oio^i 00000 .,7o/u .^.ooooo l.‘^oooo 


vvtipc n 

i VI ‘li f n 

S V Ur f i T 

VVUpf 1 ) 

cnpn2f 1 T 

SPAN! M ) 

.POO 

«1 .25/ 

,500 

1 .8/5 

1,257 

1,000 

t ‘ 500 

-,88'l 

1,500 

1 ,825 

,884 

1,000 


xj^urrli Ysuren /surrn ssurfli Vr,urfn 

y'Jsso ?, 00000 ..ottRi"? 2,0000a .b^soo 

r'lO^OS 2,00000 2.00000 1.12SO0 

r‘fi»‘i21 2,00000 70/1 1 2.00000 1 .SOOOO 


YS(l«{t) /StJiifl) .SbllHfll VSUKfl) COKUlfl) Sf’AM2(l) 

2.000 -,2o'i 2,000 ,87S 2,000 ,707 

2'*000 -.hi** 2.000 l.i7S 2,000 .S30 


FO«r.f SfcMhING LOCATltlNS TN Ni-OlREMinN 

T J »51 (j,J) YSt f l,J) VSl n,J5 

1 1 '.SOOOOOO -.IbUb/Oh 

? 1 >00»1S5 l.EiOOOOOO -,U0«H5« 

12 “'J501707 .S'OOOOOO -.7755980 

2 2 r'0827S72 T.SOOOOOO -.5^20272 


FIGURE 7.4-5. Output for POTGEM Test Case 4 (Cont’d). 
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I 3 J'«»agl35 >000000 .t 

;> 3 1>142!5‘5 1>000000 -,8e3883S 


FOKCF SkNjjI^G lOCiTlflNS IN NS-DIStr T lOfJ 


I J 

'Spf T* J) 

VS?(T, J) 

I 1 

'.5303501 

-,0000000 

^ t 

’ ISISOO*? 

r.OOOOOOO 

3 1 

'9722718 

a'.ooooooo 

1 2 

J'>37«38« 

-.0000000 

2 2 

1 ' 26lt>3t0 

r.OOOOOOO 

3 2 

+ 8TnP 
ST(iP 7/7 

r 3a582‘i,2 

2’. OOOOOOO 


,JI 

-.^ijOJSOl 
,, 3077U76 

JttN 

^,92«077h 
.,hlR7 J8U 


riGDRE 7.4-5. 


Output for POTGEM Test Case 4 (Concluded) . 
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FIGURE 7.4-6(a). Top View of FOTGEM Test Case 4. 
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FIGURE 7.4-6(b). Front View of POTGEM Test Case 4. 

Figures-76 


FIGURE 7.4-6(c). Side Viev-7 of POTGEM Test Case 4. 
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t. 

2 . 

5. 

5. 

(y, 

7. 

8. 

10 . 

n. 

12 . 

li. 

i«. 

IS. 

!<»■. 

17. 

18. 
!*». 
20 . 
21. 
22 . 
?i, 

24. 

25. 

26 . 

27, 

28. 
29. 
^0. 
?1. 

32. 

33. 

34. 

35. 

36. 


T 

T1T|.F 

TfST CASE S « SPHERE WITH S THE ClPCUMFERENTUt VARIABLE 

sktT 

•IHCBVI |C*1,COPT»0 send 
S°T1 

tINCRVJ ic»«,copt*i send 
SHTl 

iINCrVI IC*S,COPtsO send 
SPII 

»INCRV 1 IC*6,C0PT«4.V4R2(n*l sENO 
SPTl 

tl^4CRVl ic*u sEnd 

nt.EGMgNTS 

= OaTa NBPS( 1)«5,N8PV( t|»lO »EN[) 

vlbc 

cDaTA send 
SI HC 

iDaTa I0PT«? itEND 
Sl 

•^OaTA IoPTSVao SEnD 

VL 

«DaTA send 

S»i- 

iDaTA 1oPTSV*1,nTABSV»1,VAr2SV{1)*40 send 

'^UL 

mDaTa VaR2SV(i J* jBo send 

PaNi 

<OaTa send 


RqSS 

■sOaTA EAXIS(1)«0,0,J,PHI«'?0 send 

Ds.Fl 

1 B 

1 9 

IP 

1 5 


FIGURE 7.5-2. Input for POTGEM Test Case 5. 
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37. 

10 


38. 

?7 


39. 

0 


90. 

FINISH 


9J. 

iDATA 

FLTi l5*3,lUt5927,3*l,LQG(12)»T, INT(10)»1 , 1 SENQ 

92. 

STORF 

93. 

tOaTA 

IDiS sEnO 

99, 

pvtnt 


95. 

tOATA 

PRINT(n»l8*T nENn 

96. 

Slop 


FIGURE 

7.5-2. 

Input for POTGEM Test Case 5 (Concluded) . 
Figures-80 



tmWMj PAGE TB 

Pip.H)Q R 


PuTfiN UEIiMf-THir program, VERSlOM 1.^ 

«a«c«B««issc«*sassssSKsB*.«sis«s«SKSss 

TIME s 0«,‘0‘'/7f) aUlUfil 
FNTER BA^fH 

+Tmp 

TEST C;a&E S SPHfRP wJ TH S THE r iHrUMhERENT I a( VARUBLF 


vSR i \ 

♦ SRii 
»SRj1 
fSRl! 

♦ SRjt 

♦ f>Sfc GHE ^ T -> 

<-VI fit 

♦ Sf 

■t-vi: 

•fSll 

+ Vii 

♦ PANI 

+ Hnss 

♦ OSFI 
+FTNISH 
+ STf)PE 

PUP HiS BEEN aPpNtO FOR WRITING ON UN I f i 

CRFATIO<^ tT F » 0B/09/7t> 07JimS9 

creation uf GEOHETRy FTLf 


TITLF £ >EhT CASE N •- SPHtRg wITH S THE C tRCuRFE«fRTI At vARlARlt 
fU)G) s fPFFTFF*FFFFTFFFPTFFF 


FIGURE 7.5-3. Output for POTGEM Test Case 5. 
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0 


0 


0 


(INT) • o h M 5 10 

(FLTI » S.IUtS'JgT 1,0000000 

.?So6fc2eSt -02 0. 

♦PRINT 


0 0 
! .0000000 
o'. 


1 1 
r.OOOOOOO 

0. 


0 

t ,0000000 


PRJNTDU! UF GfOMFTRY FILE OaTa 


title » T^sr CASE S •« SPHERE WITH S THE CIHCUHFERENTIaL VARIABLE 


creation Tt^-f s oe/09/70 OTjjOiSO 
(iroHM) I nniitin 
(ID) * H 

(LOg) s rrrfTFFFFFFTFFFFTffF 


(IN)) * 

n 6 n 

^ lO 0 

0 0 0 

1 1 0 

n 

(f-l T ) » 

3,1911s927 

r.OOOOOOO 

1.0000000 

r.OOOOOOO 

1. 0000000 


?S'. 0662 ei 0 

' 002 S 066 

o'. 

O'. 

0 . 


0 . 


O'. 


O’. 2 S, 06628 A 


0. 


p<iNfi coRNr^ **l>ints 


I 

J 


Yf If JI 

7fT»J) 

SM f J) 

'/(T,J) 

1 

1 

-’9969173 

0000000 

,0789591 

0, 

a, 0999997 

2 

1 

-*9969173 

-.0292952 

*0796190 

16.0000000 

a ) 099999 7 

3 

1 

-.9969)73 

-.0961171 

0639797 

36>000000 

0* U9999«J7 

a 

J 

-.9969173 

-.0639797 

■0961171 

53.9999990 

0>99999 7 

R 

1 

-.99691 ’73 

-.0796190 

'.0292952 

71.9999990 

0>999997 

6 

1 

-.9969173 

-.0789591 

0000000 

90>000000 

<1 0999997 

1 

2 

-.925879S 

0000000 

.3826839 

0. 

22^5000000 

2 

2 

-.9236795 

-. 1182557 

.5639536 

IB'.OOOOOOO 

22.5000000 

3 

? 

-.9238795 

-.2299357 

'.3095979 

56>000000 

22,5000000 

u 

2 

-,9236795 

-.3095979 

.2299357 

53'.9999990 

22.5000000 

G 

? 

-.9236795 

-.3639536 

. 1182557 

71,9999990 

22.5000000 

6 

?- 

-.9236795 

-.3826839 

-.0000000 

PO’.OOOOOOO 

22.5000000 

1 

i 

-.7609060 

-.0000000 

.6999980 

o’. 

00.0999990 

? 

3 

-.7609060 

-.2006905 

>176618 

la'.ooooooo 

00|«9999»il0 

3 

3 

-, 7609060 

-.3817360 

>259195 

36'.0000000 

00,0999990 

<4 

3 

,,7609060 

-,5?S91«5 

.381 7360 

53’. 9999990 

00'* 0999990 

5 

3 

-.7609060 

-.6176618 

'.2006905 

7 r. 9999990 

00*0999990 

b 

S 

-.7609060 

-.6999980 

0000000 

Po’iooooooo 

90>999990 

1 

n 

-.5229966 

-.0000000 

.6526901 

o'. 

58 1 0999990 


FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont’d). 
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? 

M m 


.,2634803 

.8109090 

18' OOOOOOO 

58,4999990 

3 

a 0 


-.5011695 

6898004 

56.0000000 

58,4999990 

U 

u 

,S??09B6 

-.6898005 

'5011693 

55’. 9999990 

58,4999990 

S 

a 

,52?<I986 

•.6109090 

.2634803 

7i[9999990 

58.UV9999Q 

6 

u 

.5??a9(J6 

-.8526401 

-jooooooo 

90[0000000 

58[«99999() 

1 

5 

,?53a«50 

-.0000000 

.9723699 

O'. 

76.4999990 

if 

S 

.2530/450 

«, 5004786 

'.9247787 

18.0000000 

76’, 4 999990 

3 

S 

.?3300S« 

-.5715447 

. 7866658 

56’.0000000 

7b[ 4999990 

li 

s 

. 253045 a 

-.7666638 

.5715447 

53.9999990 

76.4999990 

S 

5 

,?33a05O 

-.9247787 

. 3004788 

?r,9999990 

76.4999990 

h 

S 

,?33a054 

-.9723699 

-.OOOOOOO 

90.0000000 

76[4999990 

1 

b 

, 07P4SO0 

'0000000 

’.9969173 

0. 

94.4999980 

? 


07BU59() 

-.3080644 

’.948124 7 

16^0000000 

94 [4999980 

3 

b 

078U590 

-.5859733 

‘8065251 

36.0000000 

94 [ 4999980 

a 

6 

07Ra«590 

-.8065230 

.5859733 

53‘.9999990 

94)4999980 

K 

b 

07fl4S90 

•,9481247 

. 3080644 

71'. 9999990 

94,4999980 

b 

6 

0780590 

«, 9969175 

-i OOOOOOO 

90.0000000 

90)4999980 

1 

7 

5626830 

' OOOOOOO 

’.9238790 

0. 

1 12)4999970 

ti 

7 

36268 30 

-.2854945 

'[ 878661 7 

18.0000000 

1 12.4999970 

3 

7 

3826830 

-.5430428 

. 747434,3 

36.0000000 

1 12,4999970 

u 

7 

3826834 

-.7474342 

.5430428 

S 3.9999990 

1 12,4999970 

s 

7 

3826850 

-.8786617 

. 2854945 

7 1 .9999990 

112,4999970 

b 

7 

3826830 

-.9238796 

-.0000000 

90,0000000 

1 12.4999970 

1 

b 

64 90 4(^0 

'0000000 

.7604060 

O'. 

1 30.5000000 

? 

8 

6494080 

-.2349784 

'.7231891 

18'.0000000 

liO.SOOOOOO 

3 

H 

6490480 

,,4469554 

.615181 4 

36.0000000 

1 30,5000000 

u 

b 

hOOOOBO 

-.6151814 

.4469555 

53[9999990 

1 30.5000000 


b 

bOOOObf; 

-,7?5irt9t 

'[2349784 

7l’.999999f> 

1 30,5000000 

b 

H 

b4904«n 

-.7604060 

-.0000000 

90.0000000 

1 50.5000000 

1 

q 

8526001 

’,0000000 

‘.5224967 

O'. 

148.5000000 

? 

9 

‘8526401 

-.1614610 

. 4969258 

le'.ooooooo 

3b‘.0000000 

148.5000000 

3 

<? 

[8526401 

-', 3071 1 70 

['J227105 

14H.5000000 

n 

q 

[8526401 

-.4227105 

. 3071170 

S3.994<}99o 

1 48.5000000 

S 

9 

.8526401 

-.4969257 

. 1614610 

71'. 9999990 

148.S000000 

6 

9 

‘8526401 

-,522498b 

OOOOOOO 

90[0000000 

148.5000000 

1 

1 0 

‘9725699 

' 0000000 

’,2534455 

0. 

166.5000000 

? 

tb 

9725699 

,.0721586 

[2220198 

18' OOOOOOO 

166,5000000 

3 

10 '^9723b<J9 

-.1372158 

.1888613 

36,0000000 

166,5000000 

a 

10 

.9723699 

-.1888615 

.1572158 

53.9999990 

166,5000000 


FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont'd). 
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% 

10 

■9723099 

-.2220198 

.fl7?138«» 

71.9999990 

1*6,%000000 

6 

10 

,9723099 

.,23390bS 

-.ooooooo 

90,0000000 

198.SOOOOOO 

1 

11 

1,0000000 

.0000000 

,0000000 

0, 

1(1,0000000 

160,0000000 

2 


1,0000000 

•,0000000 

.ooooooo 

160.0000000 

3 


l,OOOOOO0 

-.0000000 

0000000 

36,0000000 

160.0000000 

» 

11 

1.0000000 

-.0000000 

'ooooooo 

53,9099990 

160.0000000 

S 

1! 

r. 0000000 

*,0000000 

.0000000 

71,9999990 

1 60.0000000 

h 

11 

1.0000000 

-.0000000 

..ooooooo 

90.0000000 

160,0000000 


UNIT VE(,T1'W& 4L0NG W*KE EtfHtNTS 

I J uViNxn^j) uvwyn»Jl 

UNAVAlLAblE 


HOUNOaPV (.ONOITJON flags 
1 J BtPi.AGflfJ) 

liNAvAILAfll F 


OOUhIRT SINGUlaHITv FLaGS 
I J Dl»HAGfI»J) 


11 IN 

2 1 19 

3 1 19 

'■<1 19 

Si : p 

1 d 19 

2 2 19 

1 ? 19 

« 2 19 

S ^ 19 

13 19 

2 3 19 

3 3 19 

« 3 19 

S 5 19 


FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont'd). 
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t ^ 
? u 
i ^ 
^ li 

s « 
1 «; 
/» b 

5 *; 

t 6 
? t) 
^ b 
« b 

'i fc 
1 7 

^ 7 

i 7 

4 7 

5 7 
1 P 
? 

'?> rt 

a ft 

i; ft 

1 q 
? 9 

S 9 
a u 
S 9 

1 Ifi 

2 10 
'i' 1 1» 

a 1(1 

JO 


19 
1 9 
1 9 
1 9 
1 <^ 
1 9 
1 a 
1 9 
1 9 
1 9 
1 9 
1 9 
19 
1 9 
1 9 
19 
10 
1 9 
19 
1 0 
1 9 

IV 
1 9 
19 
19 
1" 
1 9 
19 
1 9 
1 9 
^7 
?7 
?; 
.>7 
?T 


SOlJHr^ oI:.f;U( aPITY ^i *r,5 
1 J S?^7l A(,( I , J) 


FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont'd). 
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UKAVAltASif^ 


boundary (-OMi>iTiON points 


I 

J 

» BC (1 n J ) 

VpCtI# Jl 

7fCf I, J) 

SpCf 1, J) 

VnCf I»Jl 

1 

1 


9723S99 

.,0385189 

.2305715 

8.9999990 

13,4999990 

2 

1 


9723b99 

.,1059620 

.2080013 

27.0000000 

13.4999990 

3 

1 


97?3fc99 

.,1650708 

,1650708 

45.0000000 

1 3,4999990 

a 

i 


9723S99 

.,2080015 

.1059820 

62 9999990 

15,4999990 

% 

t 


97P3(!,<»9 

..2305712 

•.0 365189 

80.9999990 

13,4499990 

1 

2 


MiiPbao? 

-.0817586 

.5160657 

8,9999990 

31 ,4999990 

2 

2 


RS2b«02 

.,2372094 

'^4655496 

27.0000000 

51.4999990 

3 

2 


«%2<*402 

.,3694^25 

.3694623 

45,0000000 

31,4999990 

u 

2 


BS26402 

.,465549b 

■.2372094 

62 9999990 

31.4999990 

s 

2 


«52640? 

-.5160657 

’.0817 368 

80^9999990 

31 ,4999990 

! 

3 


b4944ftl 

.,1189537 

.7510441 

8,9999990 

49.4999980 

? 

3 


64944«1 

.,5452171 

,8775267 

27.0000000 

49.4999980 

3 

3 


b 494 4ft 1 

.,5378882 

.5576882 

45’ 0000000 

49 4999980 

« 

3 


b494«ft! 

..8775266 

’.5452171 

62,9999990 

49%999980 

5 

5 


b49448j 

.,7510441 

.1189537 

80,9999990 

49|4999980 

t 

u 


3826H35 

.,1445266 

.9125050 

8.9999990 

67,4999980 

2 

4 


5b2t>a3S 

..4194525 

,8251827 

27' 0000000 

b7'|4999980 

3 

4 


362b«3S 

-,65328t5 

.6532815 

45.0000000 

67.4999980 

u 

4 


582toft3S 

-.8231826 

’.4194 525 

62’ 9999990 

67,4999980 


4 


582883% 

.,9125050 

’,1445268 

80 9999990 

67,4999980 

1 

5 


07H4S91 

1559522 

.9846436 

8'9999990 

85.4999980 

2 

b 


0784591 

-.4525910 

.8882598 

27.0000000 

85,4999980 

i 

h 


0784591 

-.7049270 

'.7049270 

45.0000000 

85,4999980 

a 

5 


0784591 

.,8882598 

,4525910 

62’ 9999990 

85|a999980 

s 

5 


0784591 

.,9846436 

.1559522 

80^9999990 

85,4999980 

1 

fc 


2334455 

.,1521122 

.9603985 

8.9999990 

103,4999970 


b 


2334453 

..4414467 

. 8663880 

27^0000000 

103.4999970 

i 

6 


2534453 

.,6875694 

.6875694 

45.0000000 

103,4999970 

<i 

6 


2334455 

.,8663879 

^4414467 

62‘ 9999990 

103,4999970 

5 

6 


2354453 

.,9603984 

.1521122 

80^9999990 

103.4999970 

1 

7 


5224985 

.,1333823 

.8421428 

8' 9999990 

121.4999980 

2 

7 


5224985 

.,3870905 

,7597080 

27^0000000 

121.4999980 

3 

7 


5224985 

.,6029077 

.6029077 

45.0000000 

12r.4999980 


7 


5224985 

..TS’TOSO 

.3870906 

62.9999990 

121,4999980 


FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont'd). 
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ORIGMAB PAGEl Ei 
OB POOK SUMitra 


s 

7 

'5?i«965 

-.8621428 

.1333823 

80,9999990 

I2l ,6999980 

1 

8 

.76090S9 

... )15961 

.6614523 

«[9999990 

1 39.5000000 

2 

8 

760U0S9 

-.29a«a55 

5786625 

27[0000000 

1 39.5000000 


8 

* 78090S9 

-.«59?292 

[659229? 

45'.0000000 

1 39.SOOOOOO 

ti 

8 

;7(,n«0S9 

-.5786625 

2968433 

62[9999990 

139.5000000 

5 

8 

'‘760O0S9 

-.66141523 

[ 1015961 

8o[999999C 

139'.5000000 

1 

9 

■ 9 ^ 3879 ^ 

-.0598669 

'. 3779721 

e[ 9999990 

157.5000000 

? 

9 

*9238793 

-.1 737567 

'. 360973S 

27[0000000 

157.5000000 

3 

9 

^9238793 

-.2705981 

.2705981 

45.0COOOOO 

157. 5000000 

/J 

9 

*9?38793 

,,3609735 

.1737367 

62.9999990 

157. 5000000 

S 

9 

■[9238795 

-.3779721 

0598649 

fto[9999990 

157.5000000 

1 

10 

[9989) 73 

-.0122737 

[0776932 

g[9999990 

1 75.5000000 

? 

t 0 

[99<i9t 7 3 

-.0356197 

■ 0699076 

27>000000 

1 75,5000000 

3 

t 0 

[99891 7 3 

-.0556 790 

.0556790 

4S[0000000 

175.5000000 

a 

1 0 

[996917 3 

-.0699076 

0356197 

62.9999990 

1 75,5000000 

% 

1 0 

;9W69 J 7 3 

-.0776932 

.0122737 

80.9999990 

1 75,5000000 

llivf T T 

(vJflHM '1 

3 AWn A9tfl5 





T 

J 

f f , J 1 

HNy 1 T » J ) 

UN/ M f 

Da f l» Jl 


1 

J 

-.'57?35?3 

- , Q 365668 

'.2306378 

.0223073 


2 

1 

^,9723318 

160 O 107 

.208071 3 

',0225073 


3 

t 

.,9723598 

-.1651156 

.1651154 

.0225073 



1 

9723^.12 

-.2080083 

.1060227 

'.0223073 


% 

1 

-,47?3523 

-.2506658 

.0365292 

'.0223073 


1 

2 

-,65?fo03H 

-, 081 7634 

[5161186 

'.0500647 


? 

2 

ftSe'iBOl 

-'2372028 

'.4656520 

,0500647 


3 

? 

-,8325799 

-.3695318 

5695318 

.0500647 


a 

? 

-.8525861 

-.4655978 

.2373092 

.0500647 


S 

? 

S5?9096 

-.5161673 

.08)5925 

.0500647 


1 

3 

-,6693522 

-.1189132 

.7511334 

.0731570 


? 

3 

-,6995682 

-,3453779 

>775213 

.0731370 


s 

3 

-, 6 a 99 i 57 

-.5377077 

>377078 

>731370 


a 

3 

-,6a9«7«2 

-.6775239 

.5651658 

'.0731370 


s 

3 

-,6i)947l7 

-,7510254 

*1189426 

'.0731 370 


t 

<1 

-,3827577 

-.1645023 

[9126777 

>891621 


? 


-.3828302 

-.6194507 

.8231052 

>891621 


3 


-.3827067 

-.6532747 

'.6532767 

[0891621 


ti 

« 

-.3827888 

-.8231625 

[4196153 

.0891621 




FIGURE 7 

5-3. Output 

for POTGEM Test 

Case 5 (Cont'd). 
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s 

4 . 

,3fl27S77 

-.9125005 

^1443600 

0891621 

1 

5 

0780985 

1558295 

.9846600 

0963/65 

2 

S 

,07861 i3 

..4525625 

■ 888260’ 

0963/65 

i 

s • 

,078S39B 

-.7049228 

‘ 7049228 

0963765 

a 

s 

,078«?U1 
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FIGURE 7.5-3. Output for POTGEM Test Case 5 (Cont’d). 
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FIGURE 7.5-3, Output for POTGEM Test Case 5 (Cont'd). 
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FIGURE 7.5-3. Output for POTGEM Test Case 5 (Concluded). 
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FIGURE 7.6-1. Input for POTGEM Test Case 6. 
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FIGURE 7.6-1. Input for POTGEM Test Case 6 (Cont'd). 
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FIGURE 7.6-2. Output for POTGEM Test Case 6 (Cont'd). 
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FIGURE 7.6-2. Output for POTGEM Test Case 6 (Cont'd). 
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-.500 

-.043 

.500 

.312 

1,000 

.038 

16 

' 275 

.500 

-.046 

-,500 

.275 

.275 

-.500 

-.046 

.500 

>75 

1 .000 

.036 

17 

' 240 

,50 0 

-.048 

-,500 

.240 

,200 

-.500 

-.048 

.500 

>40 

1 ,000 

,034 

> 

'2ii7 

.5 00 

-.049 

-.50 0 

.20 7 

.20 7 

-.500 

-.049 

,500 

.207 

1 .000 

,0 32 

19 

' 1 76 

,500 

-.050 

-.500 

.176 

.1/6 

-,900 

-.050 

,500 

.176 

1 ,000 

.030 

20 

'147 

.500 

-.050 

-.500 

.147 

.147 

-.500 

-.050 

.500 

>4/ 

1,000 

,027 

21 

'121 

.500 

-.049 

-,500 

.121 

.121 

-.500 

-.049 

,500 


1,000 

,025 

22 

' 0‘^7 

.500 

-.047 

-.500 

.097 

,097 

-.500 

-,047 

,500 

.097 

) ,000 

.023 

23 

' 076 

,500 

-.044 

-.500 

.076 

,076 

-,500 

-,044 

,500 

.076 

1,000 

.020 

24 

'05 7 

.500 

-.040 

-.500 

.057 

.057 

-.500 

-.040 

,500 

.057 

1 .000 

.017 




FIGURE 7.6-2 

. Output 

for POTGEM Test Case 6 (Cont 

’d). 
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as 

' ©45 

,500 

m , 036 

-.500 

.041 

,041 

-.500 

— .036 

,500 

.041 

1,000 

.015 

26 

• 0 28 

,500 

-.030 

-.500 

.028 

,028 

-.500 

.,030 

,500 

.028 

1,000 

.012 

27 

'017 

,50 0 

.,024 

— .500 

.017 

.017 

-.500 

.,024 

,500 

;017 

1,000 

.010 

a « 

’ OtiR 

,500 

.,018 

-.500 

.009 

.009 

-.500 

-,018 

.500 

’.009 

1 ,000 

.007 

S '? 

'00 5 

,500 

.,010 

.,500 

.003 

,003 

-,500 

-,oio 

.500 

.00 4 

1.000 

. 004 

50 

' OUO 

,500 

.,ooi 

-.500 

.000 

,000 

-.500 

-.003 

.500 

.000 

1.000 

.001 

51 

' 000 

■.500 

, 00fe 

-.500 

,000 

,000 

.,500 

,006 

.500 

.000 

1,000 

.001 

52 

' 0 0 5 

.500 

.015 

-.500 

.003 

.003 

-.500 

,015 

,500 

.003 

1.000 

.004 

55 

'‘ 0<'9 

.500 

.023 

-.500 

.009 

.009 

-.500 

.023 

.500 

.009 

1,000 

,007 

iu 

'017 

.500 

, 0S2 

-,500 

.017 

,017 

-.500 

.032 

,500 

.017 

1 ,000 

.010 

55 

'028 

.500 

,041 

-,500 

.028 

.028 

-,500 

,041 

.500 

.028 

1.000 

,012 

3fc 

' 04 1 

'.500 

.049 

-.500 

.04 1 

,041 

-,500 

,049 

,500 

.041 

1.000 

.015 

57 

'057 

,500 

.057 

-.500 

.057 

,057 

-.500 

.057 

.500 

.057 

1,000 

.on 

3H 

0 tb 

.500 

,065 

-,500 

, 07b 

, 07b 

-.500 

.065 

.500 

, 07b 

1,000 

,020 

5‘J 

' 0*^7 

.500 

,0 7i 

-.500 

.097 

,097 

.,500 

.073 

.500 

.097 

1,000 

,023 

40 

' 121 

.500 

,080 

»,500 

.121 

.121 

-.500 

.080 

.500 

.121 

1,000 

.025 

41 

' 147 

, 50.7 

,086 

-.500 

.147 

.147 

-,500 

, 08b 

.500 

.147 

1.000 

.027 

4 2 

' I7fc 

.500 

,091 

-.500 

.176 

, 1 7b 

-.500 

.091 

.500 

.176 

1,000 

.030 

45 

' ?0 7 

.500 

.096 

-,500 

.207 

,207 

-,500 

,096 

.500 

.207 

1,000 

.032 

44 

' 24 0 

.500 

,099 

-.500 

.240 

,240 

.,500 

.099 

.500 

.240 

1,000 

,034 

US 

'?75 

.500 

.102 

-.500 

.275 

,? 7S 

-,500 

. 10 ? 

.500 

.275 

1,000 

.0 36 

tth ■ 

'512 

'.500 

,103 

— , 5 0 0 

.312 

.312 

-.500 

.103 

.500 

.312 

1,000 

.038 

47 

' 551 

.500 

.103 

-.500 

.351 

.’■>1 

-,soo 

.101 

.500 

/351 

1.000 

,040 

48 

' 5 ' J1 

.500 

.102 

-.500 

.391 

.591 

-.500 

.102 

.500 

.591 

1,000 

.042 

4R 

' 4 34 

.500 

,100 

-.500 

.434 

,4 34 

-.500 

.100 

.500 

.434 

1,000 

, 0«3 

50 

'478 

.500 

,097 

-.500 

.478 

,478 

-.500 

.097 

.500 

478 

1 ,000 

.045 

5t 

' 5r-5 

.500 

.092 

-.500 

.523 

,523 

- ,500 

, 09 ? 

.500 

|523 

1 .000 

,046 

52 

' 15^9 ^ — 

-- . 50 - 0 — 

.086 

-.500 

.569 

.569 

-.500 

,086 

.500 

.569 

1.000 

,047 

55 

'617 

.500 

.079 

-.500 

.617 

.617 

-.500 

.079 

.500 

,617 

1 , 000 

.048 

54 

' be b 

,500 

. 07l 

-.500 

, bbb 

,666 

-, 5 q 0 

,071 

,500 

666 

1.000 

.046 

55 

' 716 

.500 

,062 

-.500 

.716 

,716 

-,500 

.062 

.500 

.716 

1,000 

,050 

5b 

' Inl 

.500 

, 05 ? 

-.500 

.767 

.767 

-,500 

.052 

.500 

.767 

1,000 

.051 

57 

'818 

.500 

.041 

-.500 

.018 

,818 

-.500 

,041 

.500 

. 8J8 

1 .000 

.051 

58 

' 8hO 

.500 

.030 

-.500 

.869 

,669 

-.500 

.030 

.500 

.869 

1.000 

. 05 ? 

5 *? 

' 922 

.500 

.018 

-.500 

.922 

.922 

-,500 

.018 

.500 

.9 22 

1,000 

.052 

bO 

'974 

.500 

.006 

-,500 

.9 74 

,974 

-, S00 

.006 

.500 

974 

1,000 

.052 


FQRct SiN-.lW(, LOCATIONS IN N } .DlRfcC T ION 

FIGURE 7.6-2. Output for POTGEM Test Case 6 (Concluded). 
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FIGURE 7.6-3(b). Oblique View of POTGEM Test Case 6. 
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t. 

a. 

3, 

4. 

s, 

<>. 

^ f 

«». 

10 . 


It, 

> 3 , 

»5. 

16, 

17 . 

18 . 
I*’. 
? 0 . 
2t. 
22. 
? 3 . 
2«. 
25. 

27 . 

20. 

29. 

30. 
51 . 

32 . 

33 , 
59. 

35. 

36. 


TtTtf 

TfST C*SE 7 ; thin MINT, WITH TWIST, CA^BER a^D DIhEONaL 
CaRY 

ThE next 2 COMMANDS Ca^SE T^E A^TS TO COjNClDE WjTh ThE TRaRING EDGE 

SRII 

.INCRVI IC« 2 , NT A8s3,VAHUl)«o., 15,479, 35,107, 

VAR2(ns-19.6025,.24,557,-3(S,411,C(JPT»l sENd 

St'Tl 

^TNCRVl IC«3,VAR2f Mao, ,3. 86975, 0.77675 «EnD 

THE NEXT command DEFINES THE TwIST *BOUT THE TWaRTNG EDGE 
As pHI * •,0097363*5**2 fnt'G.) 

St* T 1 

^tlNCRVl 1C«0 C0 PT*i.?,PARAM( 2)*0,,NTAB»3,VARI (1 )*0.,0,,-,0097363 |ENr> 

THE NEXT 2 commands Ca'JSE THE TWIST AXIS TO BE THE X.aXIS 

SRI, 

^INCPVl jr.a5,C0PTal,NTA»al,VAP2n)al, .jENtl 
T 1 

MNCRVi ICS6,COPT*0 SEND 

the next command Causes vpscaL to eouaL the local chord'. 

S: t! 

,INCRV1 IC*B,COPT«l,NTAB»3,VARlfl)»0,,l5,il79,35.107, 

VAP2( j)«j9.6o 25, lo,6l9Oj9,O,0 SENo 
T^^E next <J commands define ThE a CROSS sections 

SRI, 

UNCRVl 1C*1 l,COPT«o,SCS*o, sFNO 
SRI, 

» I NCR VI iC*l2,SCS*15,a79,CoPTs-2,PARAH{2)«0,,NTABe3, 

VA«,(l)aO,, 10, 6|9ai9, -10, 619419 ,iEND 

THE aEfine transformation Capability will re used in the npxt 
command to avoid Hand multiplying the 7/c values by ,a*c 
and to AVOjD transforming ThE independent VaRiASiE To V BY hand. 
AFTRANra),93.55896sl2*C 

S( I , 

, InCRVI IC*l3,SCSa2S,O,COPT*6,PARAM80,,NTABal5, 

VAR,(j)*-,o2»0.».0?».04i*06»,o8»,lo»«l2iil4,,,5,,,6,,16,,20,,2025,l,, 


FIGURE 7.7-2. Input for POTGEM Test Case 7. 


Figures-116 




« 0 '. 

« 2 '. 

“3. 


« 6 . 

«7. 

as. 

U9. 

so. 

52, 

53. 
5«. 

55. 

56. 

57. 

58. 
S'*. 
60 . 
61, 
62. 
63. 
bU, 
65. 
6 6 • 

67. 

68 . 
6*7. 

70. 

71. 

72. 


VAR?{ 1 )*-, 0067692,0,,. 0056767,, 00'»78fi6,, 01 30652,, 0159283. ,0170115, 
,01 79605,, 01 83928,,0183869,,0 I 83963,,0 180785,, 01 76670, 
,0176119,0., 

AFFtNF*T,iFTRANa.l.,0,,0, ,93, 55096, 1,,0, ntND 

S“T1 

ilMCRVi ICsi6^SC8*35, 107 Af TRAN(a)*72, PCFINLbT s€Nn 

? SPiNWIst segments will hE used! there will be 7 PANELS IN THf 
inboaro segment and b Panels in the outboaro segment' thert mu 

PE 5 CHORDwISE PftNELS SET BaCK 1/6 Pa.NEL FROM LEaOInS EDGE'. 

n.sE(;MENTS 

^DaTA NSEGST«?,NHP8«7,9,NBPV»S fEND 

vLrc 

iDaTA SEnO 
SLHC 

iDaTA I0PT»? 3ENd 
Sl 

“Data ioPTSv»n 5 End 

su 

»DATA lnPTSV*l,MTABSV«l,VARPSV{ nal5,679 fE ^40 

THE VLfS) rURVE OEEINEO rY THE FOLLOWING COmMaND WILL BE 
Satisfactory for the ?no segment also'. 

VL 

“Data vah2Sv(1)si, »end 

va 

tDaTA InPTSV»0 SEND 
PANl 

«DaTA SFNO 
SrGMFNT 

^DATA NsEGs»2 JEND 
VLHC 

^Oata lnPi«i sEnd 

su 

fOATA 10PTSV«1 , VaR2,«5Vo 35’ io7 nENf) 

P^Ni 

iDaTA send 

the following command eliminates the null row of Panels 


FIGURE 7.7-2. Input for POTGEM Test Case 7 (Cont'd). 
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73. 

8ETwFEN the !ST and ?N0 SEGMENTS 

7«. 

MHIg 

75. 

-.Data Ii»8 «Eno 

76. 

9pSS 

77. 

tOATA fcAXl5(l)»0,,0,,1.,PHl*«»0. sENn 

78. 

DPFL 

79'. 

1 .1 

80', 

I .1 

81 ! 

? 

82. 

1 16 

83. 

5 5 

80. 

6 

85. 

16 16 

86. 

tt 4 

87. 

a 

88. 

16 16 

89. 

5 S 

90'. 

12 

91. 

0 

92. 

FlHlsH 

93. 

^RaTA FtT(5)«l ,,0,. ,176326981 «ENO 

90. 

SlORf 

95. 

tDAtA 10«7 hENO 

98. 

PRINT 

97. 

f@ATA send 

98. 

STOP 

FIGURE 7.7-2. 

Input for POTGEM Test Case 7 (Concluded) . 
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PUT» AN ^EUM^ TPY PKnttRA’^. VFWiJiON l .i 


TIME # 0fl/0«/7A> 07 I10J I? 

enter Ha TlX 
♦TTTlf 

TEST tAi'E 7 * thin wlN(, WITH TwlSi', CaH),Ew and OIHEOPaL 


♦ cany 

♦ rHl rtT ? COMNaNdS CaMSE the axis to coincide with the trailing edge 

♦ S«I1 
+ 5R II 

♦ 1nf ..fAl tOHMaMi dEKInES THE TvilST aBOUT THE TRaJLING EPGF 

f AS F--1 * fDEG.I 

♦ SR It 

♦ tH! A UiHhaNNS Cause the TwjSi 4 X 1 S TO BE IhE X*AXivS 

♦ Sf-n 

+ SHJ 1 

♦ IHt .FiT compand CAUSES ypSC.al ^0 EOUal the LOCai lHDRI), 

♦SHU 

f IHt -Fat ii fOMHANtib define THF. .) CROSS SECTIONS 

♦ SR-Ii 

+5Rn 

■t- 1H[ --FttME THanSE OKmaIiON CAPAHjIiTY WJLL p£ USED jm T HF NEyT 

♦ tnMHfM!' TO AVOID hand MdtTlPLYJW; IHfc ;/C VALUES H¥ 1<?*L 

+ and '0 AVOJO THANSFORMrNG FHE 1 NDEPEnOEN I VAWlAHib TM V EiV hand’. 


♦ •tiTf<AN(U)a95,SSS96s:l<?*r 

fSRjl 
+ Sh 11 




c V AF. t f * ♦ t (. ) I V APp f A , j u ) T ^ 

i.OPoOiiOO -,a«73»2ii 

1.0000000 

0. 

.9800000 

’.'^9i)0'i00 ,990550a 

’.9200000 

r. 1 I0ft376 

,9000000 



I ,^^ua^8o 


',9800000 ,70«tyT9^ 

'.aaooooo 


FIGURE 7.7-3. Output for POTGEM Test Case 7. 
Figures-3.19 



,t»600o00 '.KSOOOOO .HOOOOOO l,320<»336 

'.8000000 1 . 2^20200 ' 7'»7'sOOO 1 ' 2<>805b» 0. 0. 

+ ? SfcGMENTt. WILL Pt USED. THEf^E- WILL BF ' 7 PA»^ELS IN THE 

♦ INbliAHl) SEfSMFNT ANU C) Panels TN the OUTHOaRO SfRBENT, THEHF will 
+ HE B r -ORDWf SE Panels SET HACK t /O PANfel. fHOM LEaI>1Ni, EDgE, 

+ OSfC.ME'^T ■ 

+ VLijf 
+-SLBr 
fS(, 

♦ SU 

f THE V) t- 8 , r*-J**'Vfc OEf iNEU iaY the HiUnwTNG COMMAND '^ILt. Hfc 

♦ SAln'FjCTOMY FUN THE 2ND SEGMENT aLBO". 

4VL 

♦ VII 

♦ PaNi 
fSEGMFNl 

♦ VLhL 
fSU 

♦ PANI 

+ IHf FOiLOWlNO CDMHiNp ifeUMINATpS =?HF NULL BOW UF PANELS 

+ HF TtvFE^ the 1ST AND VNR SEGMENTS 

♦ NWJ7> ■ 

t«nss 

♦ OSFl 

fHnish 

♦ S T ( ‘ H E 

FlLf 7,GM.p.-.C/LlH? Has been OPfNEn FDP WRITING ON UNIT j 

CRFAtlUN tI^E 3 0«/0«»/7P n7J5H|?rs 

creation of Geometry fue 


■.8POOOOO 1.3016S20 


TITLE s TEST case 7 . TmIN WING NITu TWIST, CAMBER AND DImEORaL 


TLOG) s 

ffFETEFFEFFFFE 

F F T F F F 






TINT) : 

0 17 6 1ft S 

0 n 0 

GOO 

0 n 




(ELT) * 

OOl’.aOORO 55, 108909 

10'.9ie9S7 

55.106999 

■.98U8077S 

0. 

1 75641818 

22U. 02290 


FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont'd). 
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, -0 1 0. o’ 0 , 


PHlwtnut UP r,pcif-ErHY fjLt DaIa 


TIUf s 

Tfsl TASt 7 - THTU WITH TWIST, CaMHER a^D 

OIXtDKAl 


CtyP AT 

II‘t 5 0«/09/76 

O7!i«l?0 




f If (iK'KJ 

= 1 1 T m n 1 





flU) * 

7 





aoto = 

‘ f F r T r T p 

F F F F F F F F 

T F F F 



fTUT) s 

" W h 

If* 9 0 

0 0 0 

0 0 0 

A 

(H n « 

■^0 1.9909000 

i>i[ 1 0^9990 

10,9189570 

55,1069990 

■ 98^8078 


??4.opp«9ort 

lo?2aoP3 

0. 

0. 

0. 

PaNU rCRr.fw pUfUTs 





! J 

V f ! , J ) 

vn,J5 

/n,4) 

Sf t , J) 

Vf f , J) 

y 1 

’9P01P49 

looooooo 

o‘ 

0. 

.9500000 

? 1 


P ' P \ J P H '-M 

.6104210 

2’. 211 269 7 

.9500000 

i 1 

a’ Hi-^'571S 

4.«2P'>7n 

1 .19597l« 

4’.4225714 

.0500000 

a 1 

fc'‘7^,9bBn 

O>3iHS70 
aifla‘51 aio 

1 .7627279 

6',6338571 

.9500000 

s 1 

H*‘h«SiaQO 

2.3167590 

b;845i4io 

.9500000 

ty 1 

1 0’ b099^aO 

1 r.O'jb'APHO 

2.H64I Mi 

1 1.0564280 

,9500000 

7 V 

IP" f.J57P.tO 

l 5lP<*7 7iao 

3'. 4 1090 12 

1 3.2677140 

.9S00000 

h ) 

ia'‘uSb5H00 

lb’ a 790000 

3'.963nul 

15.4/90000 

.9500000 

^ 1 

if*; 3 79SPB0 

171i>598880 

4.5485559 

17.6598890 

, 9500000 

1« 5 

)M;?907aoo 

1918407770 

5'. 1i«U998 

1918407770 

.9500000 

1 1 1 

pn' pib'jAbO 

ZP’ 0P16fef*0 

'* . 7257971 

2210216660 

,9500000 

IP 1 

PP" la^pqati 

24^2025S50 

6.5192591 

2412025550 

,9500000 

li 1 

0S10510 

Pb’ 3834940 

6 806/142 

2613834440 

.9500000 

1 

,-S'9BU?O0 

?«'S643i?0 

7'239i864 

28.5643310 

.9500000 

1*^ 1 

?T 91 PS900 

30l 7-)SP2V0 

7,6672207 

50.7452220 

,9500000 

lU 1 

p9’ Basoiso 

3?1 92611 00 

8.1537990 

52.92el lOG 

,9500000 

17 1 

Ml 77779P0 

3bl 1069990 

8.641 8920 

55. 1070000 

,9500000 

t ? 

alooOfa?«9 

' 0000000 

0, 

0. 

, 7500000 

? ? 

h' 570*>S19 

2,21 I28S7 

.8258557 

2l 21 128s 7 

,7500000 

i ? 

«'?395ppo 

4.422'i7l3 

1.6520706 

4.42257 14 

,7500000 

‘J P. 

9'9o5J500 

6l631Hb70 

2.423435? 

616338571 

,7500000 


FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont'd) 
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5 

2 

n' 5fe7<l920 

8' 645 1930 

3.2097096 

8‘.B951930 

. 7500000 

6 

2 

I3'22a5090 

ir 0569280 

3,9806230 

ir,0569280 

,7500000 

7 

2 

l<4'8751?b0 

15. 2677190 

9.7558605 

13.2077190 

,7500000 

H 

2 

lfe[si7»«90 

is: 9790000 

5.5350512 

15.9790000 

, 7500000 

9 

2 

1 b’ 50«8«30 

1 7 6598680 

5'.9061 1/42 

17.6598890 

.7S00000 

10 

2 

20^0926290 

19.8/407770 

6.2775055 

!9. 8907770 

.7500000 

1 1 

2 

21' Bftl 7890 

22' 0216660 

6’. 65 15689 

22.0216660 

,7500000 

12 

2 

23' 6728010 

29^2025550 

7.0306663 

29.2025550 

,7500000 

13 

2 

2S'o6O3S70 

26.3839990 

7]9325506 

26.3839990 

,7500000 

1« 

2 

?7' ?56‘>'580 

28'. 5693320 

7’,85299o8 

28.5693350 

,7500000 

I** 

2 

?9'*oi>03230 

30: 7952210 

8*2896900 

30. 7952220 

',7500000 


2 

^O* B«S<(920 

32.9261 100 

8^ 7313180 

52.9261100 

, /5OOOOO 

17 

2 

32>«17270 

35' 1069990 

9.1997380 

55. 1070000 

. 7500000 

1 

5 

8' 8?! 1?50 

'ooooooo 

0. 

0. 

.5500000 

2 

3 

10^2300160 

2,21 12857 

■^9199216 

2.2112857 

.5500000 


3 

1 r 6058790 

9;9225715 

r'‘ft25a9l2 

9'. 92257 19 

,5500000 

a 

3 

1 3'‘0537iil0 

6,6558570 

2*7200587 

6.6338571 

.5500000 

S 

3 

to'a5fe258i) 

8'^84451930 

3*60729 39 

8'.895!430 

,5500000 

h 

3 

IS’ 8S 16900 

n ‘0569280 

9*990921 7 

11,0569280 

,5'iOOOOO 

7 

5 

1 7^2582690 

13,2677190 

5*3/29569 

13.2677190 

.5500000 

H 

3 

t8’>i0l770 

15.9790000 

6*2581389 

IS',9790000 

.5500000 

9 

5 

r'O' 2756390 

1 7' 6598880 

6*9077737 

17 ’.6598890 

.5500000 

in 

3 

rl'9U 32760 

19'890777o 

6*5577993 

19,8907770 

.5500000 

11 

3 

23' 6187890 

22' 0216660 

6] 710101 5 

22.0216660 

,5500000 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Coiit’d) 
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FIGURE 7-7-3. Output for POTGEM Test Case 7 (Cont'd). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont’d) 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont'd). 
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FIGURE 7.7-3. Output for FOTGEM Test Case 7 (Cont'd). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont'd). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont’d). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont’d). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Cont’d). 
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FIGURE 7.7-3. Output for POTGEM Test Case 7 (Concluded). 
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FIGURE 7,7-4(a). Top View of POTGEM Test Case 7. 
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FIGURE 7.7-4(c). Side View of POTGEM Test Case 7. 
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